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ABSTRACT 


The  purpose  c£  the  research  reported  here  was  to 
determine  the  sources  of  scattering  of  artificial  sea 
ice,  to  compare  backscatter  measurements  semi¬ 
quant  itatively  with  theoretical  predictions,  and  to 
develop  inexpensive  polarimetric  radars  for  sea  ice 
backscatter  studies, 

In  Chapter  2,  a  brief  review  of  the  dielectric 
properties  of  sea  ice  and  of  commonly  used  surface  and 
volume  scattering  theories  is  presented^ 

In  Chapter  3,^6  provide  a  description  of  the 
backscatter  measurements  performed  and  experimental 
techniques  used- 
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^Disscussion  of  the  development  of  inexpensive  ^  ~ 
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short-range  polarimetric  radars  is  the  topic  of  g 

Chapter  4j  The  steps  taken  to  add  polarimetric 

capability  to  a  simple  FM-CW  radar  are  considered  as  ^ 

r— I 

are  sample  polarimetric  phase  measurements  of  the 
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radar.  We  also  consider  limitations  of  this  system  and  <^3 
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discuss  the  development  of  a  radar  that  avoids  the 
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oscillator  instability  problems  of  continuously  - - — - 

tuneable  oscillators  by  doing  measurements  at  discrete  C. 
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frequencies,  using  a  synthesized  frequency  source.  We 
present  results  of  simple  test  of  copolar  phase  fu*  OeJ 

capability  and  a  simple  method  of  doing  phase  ,  - 


'  calibrations  in  the  laboratory. 
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In  Chapter  5,  ve  discuss  ice  surface 

characterization  data  and  techniques,  including 
computation  of  surface  rms  height  and  correlation 
length  and  air  bubble  distribution  statistics^  We  also 


present  a  method  of  estimating  the  standard  deviation 
of  rms  height  and  correlation  length  for  cases  of  few 
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data  points. 
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:hapter  6  consists  of  comparisons  of  backscatter 
measurements  and  theory^  We  determined  that 
backscatter  from  an  extremely  smooth  saline  ice  surface 
at  C  band  ,  i.e.,  rms  height  <  0.05  cm,  cannot  be 
attributed  only  to  surface  scatter.  We  found  that  snow 
cover  had  a  significant  influence  on  backscatter  from 
extremely  smooth  saline  ice  at  C  band.  We  determined 
that  backscatter  at  C  band  from  moderately  smooth 
desalinated  ice,  at  incidence  angles  of  less  than  30 
(rms  height  ~  0.1  cm),  and  rough  saline  ice  (rms  height 
=0.5  cm)  was  dominated  by  surface  scatter.  We  found 
that  volume  scatter  from  desalinated  ice  is  important 
at  X  band  and  that  it  is  dominant  at  Ku  band  at  angles 
of  10e  or  more  from  nadir 

Conclusions  and  sugge:  >ns  for  future  research 

are  presented  in  Chapter  7. 
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Figure  2.1:  Phase  relations  for  "standard"  sea  Vce. 

Circles  on  the  brine-salt  line  indicates 
temperatures  at  which  solid  salts 
precipitate  (Assur,  1958). 

Figure  2.2:  Simplified  geometry  of  first-year  sea  ice 
( Vant  et  al . ,  1978 ) . 

Figure  2.3:  Variation  of  salinity  with  time  for  2.5  cm 
sections  of  ice  at  a  depth  between  40  and 
42.5  cm.  Broken  line  indicates  "stable" 
salinity  or  average  of  all  the  salinities 
except  the  intial  high  value  at  this  level 
(Makawo  and  Sinha,  1981). 

Figure  2.4:  Salinity  and  temperature  profiles  of  a 
typical  frazil  ice  core  (Martin,  1979). 

Figure  2.5:  Conceptual  view  of  a  brine  drainage 
channel  (Vant,  1976). 

Figure  2.6:  Temperature  and  salinity  profiles  for 
typical  multiyear  sea  ice  (McNeill  and 
Hoekstra,  1973). 

Figure  2.7:  Real  and  imaginary  parts  of  the  relative 

dielectric  constant  of  brine  (Vant  et  al., 
1978) . 

Figure  3.1:  Setup  for  CRREL’87  experiment. 

Figure  3.2:  Setup  for  CRREL'88  experiment.  Tent  (T) 

has  been  rolled  off  to  expose  pond  (P)  for 
measurements  from  gantry  (G). 

Figure  3.3:  Setup  for  CRREL'89  experiment.  University 
of  Kansas  and  Environmental  Research 
Institute  of  Michigan  radars  mounted  on 
the  gantry  (G),  University  of  Washington 
(UW)  radiometers,  and  University  of 
Massachusetts  (UM)  radar  and  radiometer 
during  measurements  on  pond  (P).  Weather 
station  (W)  is  in  foreground. 

Figure  3.4:  Vertical  growth  history  of  saline  ice 
sheet  grown  during  CRRELf87  as  seen  in 
thin  section  photographed  between  crossed 
polarizers  (courtesy  of  CRREL ) . 
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3.5:  Salinity  profiles  of  saline  ice  sheet 
grown  during  CRREL'87  experiment. 

3.6:  Growth  history  of  saline  Ice  sheet  of 

CRREL'88  experiment  as  revealed  in  thin 
sections  of  ice  structure  photographed 
between  crossed  polarizers  (courtesy  of 
CRREL). 
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weeks  of  desalination  and  growth. 
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depth  of  6  cm  in  6.5-cm  thick  ice  sheet  of 
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platelet  structure  of  crystals  at  depth  of 
5  cm  In  14.5-cm  thick  ice  sheet  of  Fig. 

3.6  after  6  days  of  growth  (courtesy  of 
CRREL ) . 

3.10:  (a)  thermally  modified  platelet  structure 
in  ice  from  depth  of  8  cm  in  19-cm  thick 
ice  sheet  of  Fig.  3.6  after  11  days  of 
growth  and  desalination  (b)  thermally 
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Fig.  3.6  after  3  weeks  of  growth  and 
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sheet  of  CRREL'88  experiment. 

i.12:  Temperature  profiles  of  CRREL'88  ice 

sheet  on  dates  measurements  were  taken. 

.13:  Estimated  dielectric  constant  profiles  at 
C  band  on  dates  measurements  were  taken 
during  the  CRREL'88  experiment:  (a)  real 
component  of  epsilon  (b)  imaginary 
component  of  epsilon. 

.14:  Two-way  absorption  loss  at  C  band  as  a 

function  of  depth  into  the  ice,  computed 
from  data  of  Fig.  3.13. 
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3.15:  Desalinated  ice  blocks  as  studied  during 
the  CRREL'89  experiment. 

3.16:  Porous  structure  of  ice  as  seen  in 

horizontal  thick  section  from  top  surface 
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3.15. 
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epsilon. 
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PAIR  (the  X-band  section  is  nearly 
identical ) . 

4.3:  Frequency-versus -dr i ving-current  response 
of  a  typical  YIG  oscillator  (Clark  and 
Swartz,  1971 ) . 
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(Clark  and  Swartz,  1971). 
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4.6: 

Effect  of  temperature  on  linearity  of  a 
typcial  heated  YIG-tuned  oscillator  (Clark 
and  Svartz,  1971 )  . 
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4.7: 

Block  diagram  of  the  Ku-band  section  of 
STEFRAD  (the  C-band  section  is  similar). 
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Vector  addition  of  returns  from  four 
targets  (Ilzuka,  1984). 
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4.9: 

Relative  phase  shift  encountered  by  r.f. 
signal  in  propagating  to  target  and  back, 
at  all  8  frequency  steps  and  various 
target  ranges  for  an  8-frequency  step 
frequency  radar. 
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:  Real  components  of  vectors  sampled  by  an 
8-frequency  step  frequency  radar  for 
various  target  ranges. 
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:  Imaginary  components  of  vectors  sampled 
by  an  8-frequency  step  frequency  radar 
for  various  target  ranges. 

Figure 

4.12 

:  Uncalibrated  copolar  phase  measured  by 
PAIR  for  several  spots  on  a  smooth  saline 
ice  surface. 

Figure 
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:  Return  from  an  8"  sphere  as  output  from 
C-band  STEFRAD  (a)  before  and  (b)  after 
removal  of  the  sy.  tem  response. 

Figure 

4 .14 

:  Return  from  an  8"  sphere  as  output  from 
Ku-band  STEFRAD  (a)  before  and  (b)  after 
removal  of  the  system  response. 

Figure 

4.15 

:  Results  of  C-band  copolar  phase 

experiment  using  inserts  of  known  phase 
shift. 

Figure 

4 . 16 

:  Results  of  Ku-band  copolar  phase 

experiment  using  inserts  of  known  phase 
shift. 

Figure 

5.1: 

Surface  profiles  for  (a)  smooth  and  (b) 
rough  saline  ice  studied  during  the 
CRREL'88  experiment. 
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Figure 

5.2: 

Surface  profiles  for  (a)  smooth  saline  ice 
sheet  and  (b)  desalinated  ice  blocks 
studied  during  the  CRREL'89  experiment. 
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3:  Computed  autocorrelations  and  minimum- 
mean-square  exponential  fits  for  surfaces 
of  (a)  smooth  saline,  (b)  rough  saline, 
and  (c)  desalinated  ice  sheets  studied 
during  the  CRREL'88  experiment. 

4:  Computed  autocorrelations  and  minimum- 
mean-square-exponential  fits  for  surfaces 
of  (a)  smooth  saline  ice  sheet,  and  (b) 
desalinated  ice  blocks  studied  during  the 
CRREL'89  experiment. 

5:  Computed  autocorrelations  and  ±.  1  standard 
deviation  limits  obtained  from  Monte  Carlo 
simulations  for  surfaces  of  (a)  smooth 
saline,  (b)  rough  saline,  and  (c) 
desalinated  ice  sheets  studied  during  the 
CRREL’88  experiment. 

6:  Computed  autocorrelations  and  ±  1  standard 
deviation  limits  obtained  from  Monte  Carlo 
simulations  for  surfaces  of  (a)  smooth 
3aline  ice  sheet,  and  (b)  desalinated  ice 
blocks  studied  during  the  CRREL’89 
exper iment . 

7:  Bubbles  in  a  3-mm  thick  section  taken  from 
the  top  of  a  sample  of  a  desalinated  ice 
block  studied  during  the  CRREL'89 
experiment,  as  visually  identified  in  a 
thick-section  photograph. 

8:  Autocorrelation  of  air  bubbles  in  (a)  a  3- 
mm  thick  horizontal  section  of  the  top  and 
(b)  3-mm  thick  vertical  section  near  the 
middle  of  a  sample  taken  from  the 
desalinated  ice  blocks  studied  during  the 
CRREL'89  experiment,  computed  from  thick- 
section  photographs  that  were  processed  as 
shown  in  Fig.  5 . 7  .  (Courtesy  of  CfiREL) 

9:  (a)  roughness  and  (b)  computed 

autocorrelation  of  the  surface  of  the  snow 
covering  the  saline  ice  sheet  on  Feb.  3, 
1989  . 

1:  Returns  from  target  and  system  as  seen  in 
range-dlscr iminated  outputs  of  (a)  X-band 
SOURCESCAT  (b)  Ku-band  STEFRAD. 


Figure  6.2:  Relative  power  return  at  X  band  from  an  8" 
Luneberg  lens  as  the  bistatic  angle,  6,  is 
varied  (courtesy  of  P.  Racette  and  R. 
Forster  (19893). 

Figure  6.3:  Backscatter  cross-section  of  smooth  (i.e., 
rms  height  <  0.1  cm)  saline  ice  sheet  at  X 
band,  as  determined  from  data  obtained 
during  the  CRREL'87  experiment. 

Figure  6.4:  (a)  measured  backscatter  from  artificially 
grown  sea  ice  at  X  band  and  30  degrees 
incidence,  contrasted  with  point  target 
return  (b)  simulation  of  relative 
contributions  of  surface  and  volume 
scatter  to  backscatter  from  artificially 
grown  sea  ice  at  X  band  and  30  degrees 
incidence . 

Figure  6.5:  Backscatter  cross-section  of  textured 

saline  ice  sheet  at  X  band,  as  determined 
from  data  obtained  during  the  CRREL'87 
experiment . 

Figure  6.6:  Effect  of  variation  in  (a)  correlation 

length  and  (b)  rms  height  on  predictions 
of  the  Small  Perturbation  Method  at  X 
band . 

Figure  6.7:  Backscatter  cross-section  of  smooth  saline 
ice  sheet  at  C  band,  as  determined  from 
data  obtained  during  the  CRREL'89 
experiment,  compared  with  Small 
Perturbation  Method  predictions. 

Figure  6.8:  Backscatter  cross-section  of  smooth  saline 
ice  sheet  at  C  band,  as  determined  from 
data  obtained  during  the  CRREL'89 
experiment,  compared  with  predictions  of  a 
simple  layer  model. 
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Figure  6.9:  Effect  of  snow  cover  as  seen  in  range- 

discriminated  output  of  X-band  SOURCESCAT 
during  the  CRREL'87  experiment.  (a) 
return  from  bare  smooth  saline  ice  surface 

-  peak  at  slant  range  of  about  3.6m 
corresponds  to  surface,  (b)  broadening 
effect  caused  by  volume  scatter  within  15- 
cm  snow  cover,  and  (c)  return  after  15  cm 
of  loose  snow  packed  to  thickness  of  7.5cm 

-  peak  at  3.4m  corresponds  to  snow  surface 
and  peak  at  3.8m  corresponds  to  ice 
surface . 

Figure  6.10:  Effect  of  6.5cm  cover  on  backscatter 

cross-section  of  smooth  saline  ice  at  X 
band,  as  determined  from  data  obtained 
during  the  CRREL'87  experiment,  for  (a) 

VV  polarization  and  (b)  HH  polarization. 

Figure  6.11:  Effect  of  7cm  snow  cover  on  backscatter 
cross-section  of  smooth  saline  ice  at  C 
band,  as  determined  from  data  obtained 
during  the  CRREL'89  experiment  (see  data 
points  on  Figs.  6.7  and  6.13). 

Figure  6.12:  Backscatter  cross-section  of  snow  covered 
smooth  saline  ice  surface  at  C  band,  as 
determined  from  data  obtained  during  the 
CRREL'89  experiment,  compared  with 
predictions  of  the  layer  model  of  Ulaby 
et  al.  (19821  . 

Figure  6.13:  Backscatter  cross-section  of  snow  covered 
smooth  saline  ice  surface  at  C  band,  as 
determined  from  data  obtained  during  the 
CRREL'89  experiment,  and  predicted 
relative  contributions  of  snow  surface, 
snow  volume,  and  ice. 

Figure  6.14:  Comparison  of  backscatter  cross-sections 
of  smooth,  rough,  and  desalinated  ice  at 
C  band,  as  determined  from  data  obtained 
during  the  CRREL'88  experiment  (see  data 
points  on  Figs.  6.15  and  6.16). 

Figure  6.15:  Backscatter  cross-section  of  rough  saline 
ice  sheet  at  C  band,  as  determined  from 
data  obtained  during  the  CRREL'88 
experiment,  compared  with  predictions  of 
the  Small  Perturbation  Method. 
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Figure  6.16:  Backscatter  cross-section  of  desalinated 
ice  sheet  at  C  band,  as  determined  from 
data  obtained  during  the  CRREL'88 
experiment,  compared  with  predictions  of 
the  Small  Perturbation  Method. 


Figure  6.17:  Effect  of  temperature  on  backscatter  from 
desalinated  ice  at  (a)  C  band  and  (b)  X 
band,  as  determined  from  data  obtained 
during  the  CRREL'88  experiment. 

Figure  6.18:  Backscatter  cross-section  of  desalinated 
ice  at  (a)  C  band  and  (b)  Ku  band,  as 
determined  from  data  obtained  during  the 
CRREL'89  experiment,  compared  with 
combined  predictions  of  Small 
Perturbation  Method  surface  and  Effective 
Medium  volume  scattering  theories. 

Figure  6.19:  Backscatter  cross-section  of  desalinated 
ice  at  (a)  C  band  and  (b)  Ku  band,  as 
determined  from  data  obtained  during  the 
CRREL'89  experiment,  compared  with 
predictions  of  the  Effective  Medium 
theory  at  three  different  values  of 
volume  fraction  of  air. 


Figure  6.20:  Backscatter  cross-section  of  desalinated 
ice  at  (a)  C  band  and  (b)  Ku  band,  as 
determined  from  data  obtained  during  the 
CRREL'89  experiment,  compared  with 
predictions  of  the  Effective  Medium 
theory  for  three  different  values  of  the 
imaginary  component  of  the  dielectric 
constant . 


Figure  6.21:  Backscatter  cross-section  of  desalinated 
ice  at  (a)  C  band  and  (b)  Ku  band,  as 
determined  from  data  obtained  during  the 
CRREL'89  experiment,  compared  with 
predictions  of  the  Effective  Medium 
theory  for  three  different  values  of 
correlation  length  of  spacing  between  air 
bubbles  . 
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1.0  Introduction 

The  Arctic  region  has  attracted  considerable 
attention  in  recent  years  for  several  reasons.  The 
study  of  sea  ice  as  an  indicator  of  global  climate  is 
among  the  most  important  scientific  reasons. 
Commercially,  the  Arctic  offers  a  wealth  of  mineral 
resources  which  will  undoubtedly  be  exploited  as  those 
in  the  temperate  regions  become  scarce. 

This  region  has  remained  largely  unexploited, 
however,  because  of  its  uninviting  environment.  Weather 
can  be.  severe,  sunlight  is  absent  for  several  months, 
and  in  some  areas  sea  ice  is  always  present,  making 
navigation  hazardous  and  often  impossible. 

A  knowledge  of  the  morphology,  extent,  and  dynamics 
of  sea  ice  is  essential  to  successful  exploration  and 
development  of  the  Arctic  and  Antarctic  regions  and  to 
better  understanding  of  long-term  global  climatic 
changes.  Since  the  long  polar  winters  and  frequent 
cloud  cover  limit  the  utility  of  visible  and  infrared 
sensors,  microwave  remote  sensing  has  long  been  used  as 
an  ice-surveillance  tool.  Microwave  systems  can  see 
through  clouds  and  darkness  and  are  thus  suited  for  this 
task  . 

Microwave  sensors  fall  into  two  major  categories: 
passive  sensors  that  monitor  inherent  black-body 
radiation;  active  sensors  that  monitor  reflected 


radiation  [Thomas,  19861.  Because  of  the  easy  access  to 
passive  microwave  data  fr  satellites,  radiometers  have 
been  well  tested  for  monitoring  areal  ice  concentration 
and  surface  wetness.  Radar,  especially  synthetic 
aperture  radar  ( SAR ) ,  is  generally  more  suited  for 
determining  ice  types,  surface  features,  and  movement, 
because  of  its  finer  resolution.  The  principal 
disadvantage  of  SAR  is  the  high  rate  of  data 
acquisition.  Radar  altimeters  may  also  provide 
estimates  of  the  ice  edge  and  areal  ice  concentration. 
Much  research  currently  underway  is  aimed  at  determining 
what  information  can  be  obtained  from  microwave  sensors, 
optimum  parameters  for  the  sensors,  and  methods  of 
automatically  interpreting  the  data  gathered  by  the 
sensors.  This  research  is  of  great  importance  in  light 
of  the  present  operation  of  the  Soviet  Kosmos  1870  SAR 
and  the  approaching  deployment  of  the  European  Space 
Agency  (ESA)  Earth  Remote  Sensing  (ERS-1)  satellite, 
slated  for  launch  in  the  early  1990's,  Canadian 
RADARSAT,  Japanese  JEP.S-1,  and  the  Earth  Observation 
System  (EOS),  scheduled  for  the  mid-1990's. 

1.1.0  Background 

1.1.1  History  of  remote  sensing  of  sea  ice 

Radar  was  used  for  experimental  mapping  of  Arctic 
ice  as  early  as  1962,  when  the  U.S.  Army  Cold  Regions 
Research  and  Engineering  Laboratory  ( CRREL )  flew  a  Side- 
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Looking  Airborne  Radar  ( SLAR )  over  Arctic  pack  ice 
[Anderson,  1966].  In  his  analysis  of  the  data,  Anderson 
shoved  that  ice  types  can  be  identified  from  SLAR 
imagery.  To  verify  the  potential  of  radar  to 
discriminate  between  ice  types,  the  National  Aeronautics 
and  Space  Administration  (NASA),  Naval  Oceanographic 
Office,  and  CRREL  conducted  a  joint  experiment  over  the 
Arctic  with  a  13.3  GHz  scatter ometer  in  1967  [Rouse, 
1969].  In  1968,  the  Naval  Oceanographic  Office  flew  an 
X- , C- , L- , P-band  mult i frequency  SLAR  over  the  ice  fields 
north  of  Alaska  [Ketchum  and  Tooma,  1973].  In  their 
report,  Ketchum  and  Tooma  concluded  that  X-band  holds 
the  greatest  potential  for  sea  ice  mapping,  determining 
the  stage  of  ice  development  and  for  fracture-  pattern 
analysis.  They  also  concluded  that  L-band  sees  only 
more  prominent  features  such  as  ridges  and  hummocks.  In 
1969,  the  U.S.  Coast  Guard  conducted  ice  mapping 
experiments  in  the  Northwest  passage  using  a  Ku-band, 
unfocussed  synthetic-  aperture  radar  (SAR)  [Johnson  and 
Farmer,  1971].  In  their  report,  Johnson  and  Farmer 
concluded  that  the  ice  concentration,  ice-floe  size  and 
number,  and  water  openings  are  readily  determined  from 
SLAR/SAR  images.  They  further  concluded  that,  through 
careful  analysis,  the  age  of  ice,  drift,  surface 
topography,  fractures,  and  pressure  ridge 
characteristics  can  be  discerned  to  varying  degrees. 
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Others  involved  in  the  early  years  in  investigating 
the  potential  for  microwave  radar  as  a  sea  ice  mapping 
tool  include  the  Soviets,  who  in  the  early  1970's, 
reported  on  sea  ice  measurements  conducted  with  the 
TOROS  SLAR  [Glushkov  and  Kamarov,  1971;  and  Loshchilov 
and  Voyevodin,  1972]. 

Operational  use  of  such  radars  started  in  the  early 
1970's.  For  example,  airborne  and  spaceborne  radars 
have  been  used  to  monitor  ice  on  the  Great  Lakes 
[Hagman,  1976;  and  Gedney  et  al.,  1975],  and  along  the 
Alaskan  coast  [Gedney  et  al.,  1977]  to  assist 
navigators;  and  the  Soviets  reportedly  used  data  from 
the  TOROS  SLAR  to  provide  sea  ice  maps  to  navigators  and 
fishermen  [Loshchilov  et  al.,  1978].  To  improve  their 
operational  capability  many  field  experiments  were  and 
are  being  conducted,  during  the  1970's  and  1980's,  in 
the  Arctic  and,  more  recently,  in  the  Antarctic.  Most 
experiments  were  conducted  along  ice  edges  or  in 
marginal  ice  zones  because  the  processes  in  these 
regions  impact  importantly  on  global  climate  and  also 
greatly  influence  oil  and  mineral  exploration,  naval 
operations  and  commercial  fishing  (Luther,  1988].  Major 
field  experiments  of  the  1970's  included  the  Arctic  Ice 
Dynamics  Joint  Experiment  (AIDJEX)  in  the  Beaufort  Sea 
(Gray  et  al.,  1977;  and  Campbell  et  al.,  1978],  the 
joint  U.S.-  U.3.S.R.  Bering  Sea  Experiment  (BESEX) 


[ Gloersen  et  al.,  1975;  and  Ramseier,  19751,  and  the 
Norwegian  Remote  Sensing  Experiment  in  the  Norwegian  Sea 
[Matzler  et  al.,  1984;  and  Svendsen  et  al.,  1983]. 

Similar  experiments  conducted  in  the  1980's  include  the 
Marginal  Ice  Zone  Experiment  (MIZEX)  in  the  Bering  and 
Greenland  Seas  ([Luther,  1988;  and  Shuchman  et  al., 
1988]),  the  Labrador  Ice  Margin  Experiment  (LIMEX)  in 
the  Labrador  Sea  [McNutt  and  Raney,  1988],  and  the 
Bothnian  Experiment  in  Preparation  for  ERS-1  (BEPERS)  in 
the  Baltic  Sea  [Hallikainen  et  al.,  1988].  The  role  of 
remote  sensing  in  these  experiments  is  two-fold:  1)  to 
provide  synoptic  data;  and  2)  to  improve  understanding 
of  interaction  between  electromagnetic  radiation  and  and 
the  ice/ocean  surface. 

Because  field  conditions  have  been  found  to  be 
highly  variable,  a  laboratory  program  was  established  at 
the  U.S.  Army  Cold  Regions  Research  and  Engineering 
Laboratory  ( CRREL )  in  Hanover,  New  Hampshire  in  the 
early  1980's  to  examine  the  interaction  between 
electromagnetic  radiation  and  sea  ice. 

Research  on  microwave  remote  sensing  of  sea  ice  has 
proceeded,  during  the  1970's  and  1980's,  along  a  number 
of  fronts:  (1)  determination  of  the  scene  parameters 
that  can  be  extracted  from  backscatter  measurements;  (2) 
development  of  procedures  for  automatically  extracting 
scene  parameters;  and  (3)  improvement  of  our 
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understanding  of  microwave  interactions  with  sea  ice. 

The  latter  is  considered  in  some  detail  in  Chapters  2 
and  6  . 

1.1. 1.1  Determination  of  scene  parameters 

Sea-ice  parameters  that  are  desirable  to  measure 
i nclude : 

1)  Thickness 

2)  Areal  concentration 

3)  Ice  type 

4)  Surface  conditions 

5)  Floe  size  and  movement 

6)  Location  of  leads  and  polynyas 

7)  Location  of  the  ice  edqe 

Radio-frequency  techniques  have  been  used  to 
determine  ice  thickness  (and,  thus,  ice  type).  Knowing 
of  the  application  of  r.f.  techniques  to  study  glaciers, 
v.  Bogorodsky  and  V.  Rudakov  suggested  in  1956  that 
similar  methods  be  used  for  measurements  of  the 
thickness  of  sea  ice  and  fresh  water  ice  (Bogorodsky, 
19781.  Finkelstein  and  Kutev  [1972]  used  video  pulses 
of  nanosecond  duration  and  repetition  rates  of  a  few 
tens  of  MHz  and  found  that  ice  with  salinity  of  10  - 
13o/oo  could  be  penetrated  to  depths  of  .75  to  1.5m. 
Wadhams  et  al.  [1987]  reported  on  the  use  of  a  100  MHz 
impulse  radar  to  determine  ice  thickness  in  the  Atlantic 
sector  of  the  Antarctic  Ocean.  VHF  techniques  are  of 
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limited  usefulness  from  spacecraft,  however,  because 

adequate  gains  are  not  achieveable  using  reasonable-size 
antennas,  and  because  of  ionospheric  propagation 
effects.  Since  GHz- frequency  signals  penetrate  only  a 
few  cm  into  first-year  ice  and  undergo  considerable 
volume  scatter  in  multiyear  ice,  they  are  not  useful  for 
direct  thickness  measurements  except  for  very  thin  ice. 

The  topic  of  areal  ice  concentration  has  been 
thoroughly  addressed.  Many  papers  have  been  written  on 
concentration  estimation  from  passive  microwave  data; 
e.g.,  [Alfultis  and  Martin,  1987];  [Carsey,  1982]; 
[Cavalieri  et  al.,  1986];  [Gloersen  and  Campbell,  1988]; 
[Comiso  and  Sullivan,  1986a];  and  [Comiso,  1986b], 

Others  were  directed  toward  concentration  estimation 
from  radar  data;  e.g.,  (Loschilov  et  al.,  1978]; 

[Comiso,  1986b];  and  (Wackerman  et  al.,  1988].  A  few 
papers  directly  compare  the  two  methods;  e.g.,  [Burns  et 
al.,  1987];  [Martin  et  al.,  1987];  and  [Steffen  and 
Maslanik,  1988],  The  accuracy  of  the  estimates  of  ice 
and  open  water  have  been  reported  to  be  between  5%  and 
13%  depending  on  whether  the  measurements  were  near  the 
ice  edge  or  in  the  ice  pack.  The  schemes  used  to  obtain 
the  estimates  consisted  of  a  ratio  of  polarization 
ratios  for  the  passive  devices,  and  modal  analysis  of 
histograms  of  local  statistics  for  the  active  devices. 

The  problem  of  classifying  ice  types  in  winter 
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conditions  using  microwave  emission  and  backscatter 
measurements  occupied  the  attention  of  many 
experimenters.  Rouse  [1969]  indicated  the  potential  use 
of  13.3  GHz  scatterometer  data  to  discriminate  between 
old  and  new  ice.  In  1971,  as  a  part  of  the  Arctic  Ice 
Dynamics  Joint  Experiment  (AIDJEX),  airborne 
measurements  of  microwave  emission  between  0.3  and  3.0cm 
were  combined  with  extensive  surface-comparison 
measurements  to  further  investigate  similar  potential 
for  passive  devices  [Gloersen  et  al.,  1973].  Using 
rather  crude  images  from  a  16.5  GHz  SLAR,  Parashar 
[1974]  classified  ice  types  according  to  image  intensity 
as:  (a)  open  water,  (b)  thin  first-year  ice,  (c)  thick 
first-year  ice,  and  (d)  multiyear  and  very  thick  first- 
year  ice.  Similarly,  using  SLAR  data,  Loschilov  et  al. 
[1978]  were  able  to  classify  ice  types  according  to  (a) 
new  ice  and  nilas,  { b )  pancake,  (c)  grey,  (d)  grey- 
white,  and  (e)  multiyear  ice.  From  both  13.3  GHz 
scatterometer  and  SLAR  data  Gray  et  al.  [1977] 
determined  that  the  typical  contrast  between  first-year 
and  multiyear  ice  is  B  -  10  dB  for  like  polarization  and 
15  -  18  dB  for  cross  polarization.  More  recently  it 

has  been  reported  that  emission  data  at  33.6  GHz  can  be 
used  to  distinguish  between  open  water,  frazil,  old  ice, 
and  young  first-year  ice  [Eppler  et  al.,  1986]. 

Using  data  collected  by  a  surface-based 
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scatterometer  in  the  Beaufort  Sea,  Onstott  et  al. 

confirmed  that  <r°(0)  can  be  used  to  discriminate  thick 
first-year  (FY)  ice  from  multiyear  (MY)  ice  [Onstott  et 
al.,  1979a;  Onstott  and  Delker,  1979b;  and  Onstott, 

1980].  They  also  confirmed  that  the  contrast  between  MY 
and  FY  ice  is  greater  for  cr oss -po lar izat i on  than  it  is 
for  like-polarization.  Kim,  Onstott  and  others 
concluded  that  the  most  effective  discrimination  of  FY 
from  MY  occurs  for:  (a)  microwave  frequencies  between  8 
and  18GHz,  and  (b)  incidence  angles  greater  than  30°. 

The  least-effective  frequencies  studied  were  between  1 
and  2GHz  [Kim,  1984a;  Kim  et  al.,  1985;  Onstott  et  al., 
1979a;  Onstott  and  Delker,  1979b;  and  Onstott,  19801. 

A  combination  of  frequencies  can  reduce  ambiguity. 
Using  400MHz  and  13.3GHz  data  collected  over  Alaska 
during  a  NASA  Earth  Resources  Aircraft  Program  mission, 
Parashar  (19741  determined  that  the  400MHz  data  helped 
to  eliminate  the  MY/thin-FY  ice  ambiguity  observed  in 
the  13.3GHz  data.  Cluster  analysis  of  brightness 
temperatures  at  two  and  three  different  frequencies  has 
also  proven  to  be  useful  in  discrimination;  e.g., 

[Comiso  and  Sullivan,  1986al;  and  [Comiso,  1986b]. 

Using  18  and  37  GHz  radiometer  data  obtained  by  the  NASA 
CV-990  over  the  Bering  Sea,  Cavelieri  et  al.  [1986] 
were  able  to  distinguish  among  new,  young,  and  first- 
year  ice.  Rothrock  et  al.  [19881  did  a  principal- 
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component  analysis  on  10  channels  of  sea  ice  data 
collected  by  the  Nimbus  7  SMMR  and  determined  that  2  or 
3  of  the  channels  explained  most  of  the  data. 

Some  experimenters  have  combined  simultaneous 
active  and  passive  microwave  measurements  in  attempts  to 
improve  discrimination  between  ice  types.  Gray, 

Hawkins,  and  Livingstone  combined  emissivity  data  at 
19.4GHz,  H-po lar i za t i on  and  45°  incidence,  with  13.3GHz 
like-  and  cross-polarized  backscatter  data  [Gray  et 
al.,  1982;  Hawkins  et  al.,  1980;  and  Livingstone  et  al., 
1980].  They  concluded  that  a  combined  scheme  improved 
classification  of  sea-ice  types  and  conditions. 

Kurskaya  and  Kutuza  [1984]  stated  that  radiometric 
measurements  aided  in  interpreting  SLAR  images. 

Experimenters  have  determined  that  ice 
classification  is  possible  during  winter  months;  the 
problem  is  much  more  difficult  during  the  summer  season 
due  to  surface  melting.  Using  data  from  the  Mould  Bay 
1982  experiment.  Nimbus  5,  and  Seasat  ( July-October 
1978),  Carsey  concluded  that  the  summer  values  of 
backscatter  and  brightness  temperature  are  controlled  by 
surface  morphology  and  wetness  [Carsey,  1985  ].  L-,C-,X-/ 

and  Ku-band  scatterometer  data  collected  at  Mould  Bay 
1982  by  Onstott  and  Gogineni  showed  reductions  and  even 
reversals  in  contrast  between  FY  and  MY  ice  because  of 
the  melting  [Onstott  and  Gogineni,  1985;  and  Gogineni, 
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1984].  They  concluded  that  lower  frequencies,  such  as 

L-band,  are  slightly  better  for  discriminating  between 
ice  types  during  the  melt  season.  In  a  nutshell,  a 
fully  reliable  method  of  classifying  ice  by  remote 
sensing  techniques  during  summer  has  not  been  found. 

Ambiguity  in  interpretation  can  also  occur  in  the 
winter  months,  but  not  to  the  same  extent.  For  example, 
thin  ice  can  be  misinterpreted  as  multiyear  ice 
[U. Kansas],  and  snow  cover  can  greatly  alter  the 
backscatter  from  sea  ice  [Kim,  1384b;  and  Lohanick  and 
Grenfell,  1986]. 

The  ability  to  track  moving  ice  is  important 
because  of  the  hazards  it  presents  and  to  map  Arctic 
ocean  currents.  Curlander  et  al.  [1985]  used  an  ice- 
feature  tracking  scheme  to  study  ice  movement  between 
successive  Seasat  SAR  images.  They  determined  mean 
displacements  of  15.3km/day  in  the  marginal  ice  zone  and 
5km/day  in  the  central  ice  pack  with  an  overall  error  of 
150-200m.  Fily  and  Rothrock  [1987]  used  a  scheme  of 
nested  correlations  to  track  sea  ice  movement.  Their 
technique  consisted  of  correlations  at  several  levels  cf 
resolution  from  coarse  to  fine.  The  fine-resolution 
correlations  are  driven  by  the  results  of  coarse- 
resolution  correlations.  Ninnis  et  al.  [1986],  and 
Collins  and  Emery  [1988]  sought  to  optimize  the 
correlation  technique  by  using  a  matched-filtering 
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approach.  Vesecky  et  al.  [19861  discussed  the  use  of 
both  the  ice-feature  tracking  technique  and  the  nested- 
correlation  technique.  Shuchman  et  al.  [1986]  used 
sequential  SAR  images  to  study  the  interactions  of  sea 
ice,  ocean,  and  atmosphere.  Fily  and  Rothrock  [1986] 
studied  lead-area  changes  in  successive  SAR  images. 

The  ability  to  delineate  geometric  shapes  can  also 

be  advantageous.  For  example,  ridges  and  icebergs  can 

be  hazardous  to  navigation.  Askne  and  Johansson  [1988] 

tried  to  determine  the  ridge  characteristics  that  can  be 

discerned  from  radar  by  using  a  stochastic  triangular 

model  for  ridges.  Larson  et  al.  [1978]  considered  the 

use  of  SAR  systems  for  icebe.g  detection  and 

characterization.  Using  data  collected  by  an  X-/L-band 

o  o 

SAR  at  depression  angles  of  ’  and  30  they  concluded 
that  it  is  possible  (1)  to  detect  icebergs  in  a  variety 
of  sea  clutter,  (2)  to  identify  icebergs  by  shape,  and 
(3)  to  estimate  size  from  shadow  information.  Gray 
[1980]  similarly  concluded  that  large  incidence  angles 
enhance  the  detectability  of  ridges  and  other  objects, 
such  as  icebergs,  in  an  open  sea  or  when  surrounded  by 
sea  ice.  Shaw  et  al.  [1988]  used  data  obtained  off  the 
coast  of  Newfoundland  to  study  for  example,  the 
relationship  of  iceberg  responses  to  incidence  angle  and 
the  look  direction  relative  to  predominant  wave 
direction.  Wadhams  [1988]  used  shipboard  radar  to 
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monitor  icebergs  in  the  south  Atlantic,  and  Walsh  et  al. 


[1986]  determined  that  ground-wave  radars  are  effective 
ice-hazard  sensors.  Marthaler  and  Heighway  [1978] 
discussed  radar  image  processing  of  SLAR  data  for 
detection  and  identification  of  icebergs. 

Tn  recent  years,  considerable  emphasis  has  been 
placed  on  polarimetric  measurements:  (1)  to  determine  if 
useful  sea-ice  target  information  is  contained  in  the 
copolar  phase;  and  (2)  so  that  targets  can  be  further 
studied  via  the  technique  of  polarization  synthesis. 
During  the  last  few  years  such  measurements  have  been 
made  at  CRREL;  e.g.,  [Onstott  et  al.,  1988].  In 
addition,  the  Jet  Propulsion  Laboratory  (JPL)  has 
recently  imaged  sea  ice  regions  north  of  Alaska  with  a 
P - , L - ,  and  C-band  polarimetric  SAR.  Analysis  of  these 
images  is  still  in  progress. 

1.1.1. 2  Automat ica 1 ly  extracting  scene  parameters 

The  routines  discussed  in  the  previous  section  for 
tracking  ice  movement  were,  out  of  necessity,  at  least 
partly  automated.  Considerable  work  has  also  been  done 
on  developing  automated  routines  for  ice-type 
classification.  Holmes  et  al.  [1982]  developed  a 
digital  procedure  for  ice-type  discrimination.  Their 
method  reduces  the  speckle  effect  frequently  found  in 
SAR  images,  without  visibly  degrading  the  edges  of  the 
different  ice  types.  It  uses  a  nonlinear  isotropic 
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filtering  technique  called  neighborhood  processing. 

Lyden  et  al.  [1982]  considered  the  importance  of 
textural  differences  between  ice  types.  Having  found 
evidence  that  the  image  statistics  of  FY  ice  differ  from 
those  of  MY  ice,  they  considered  hypothesis-testing 
techniques  for  the  different  probability  distributions 
of  different  ice-types.  Fily  et  al.  [1986]  discovered 
that  they  could  reliably  discriminate  open  water  from 
thin  ice  via  a  bivariate  distribution  of  brightness  and 
variance.  Wackerman  et  al.  [1988]  presented  routines  to 
determine  ice  concentration  and  to  classify  sea-ice 
types.  The  ice-  concentration  algorithm  is  an  iterative 
scheme  that  fits  the  observed  histogram  to  the  weighted 
sum  of  the  probability  density  functions  (pdfs)  of  open 
water  and  ice.  The  ice-  class i f icat i on  routine  consists 
of  a  modal  analysis  of  histograms  of  local  statistics. 
Bogorodsky  and  Darovskih  [1984]  and  Murthy  and  Haykin 
[1987]  developed  Bayesian  classification  schemes. 

Murthy  also  discussed  nonparametr ic  classification  and 
multilevel  thresholding  techniques.  Parashar  [1974] 
tested  a  similar  technique,  developing  a  decision  rule 
on  the  basis  of  information  in  a  training-data  set. 

Using  7  categories,  he  found  that  this  routine 
interpreted  66%  of  the  remaining  data  set  correctly. 
Raizer  et  al.  [1986]  used  numerical  models  of  emission 
characteristics  to  retrieve  thickness,  temperature. 
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salinity,  and  density.  The  program  "iceman"  has  been 

developed  at  the  Jet  Propulsion  Laboratory  to  classify 
ice  types,  and  to  study  ice  surface  texture  and  ice 
motion  [Lee,  1988], 

1.2  Organization 

This  thesis  is  organized  into  7  chapters.  A  brief 
outline  of  each  is  presented  here. 

An  overview  of  sea  ice  structural  and  dielectric 
properties  and  a  review  of  scattering  theories  being 
used  to  interpret  microwave  backscatter  from  sea  ice  are 
the  subjects  of  Chapter  2. 

Chapter  3  consists  of  a  description  of  backscatter 
measurements  and  surface  characteristics  obtained  during 
three  years  of  experiments  at  CRREL .  Backscatter 
measurement  techniques  are  described,  including  spatial 
sampling  and  calibration  procedures.  Emphasis  is  placed 
on  surface  characterization,  including  a  description  of 
temperature  and  salinity  profiles,  surface  roughness 
parameters,  inhomogeneities,  and  snow  cover  parameters. 

Radar  system  developement  and  performance  is  the 
subject  of  Chapter  4.  FM-CW  radars  are  commonly  used  as 
short-range  scatter ome ter s .  Even  with  high  quality  YIG 
oscillators,  the  performance  of  these  radars  is  limited 
due  to  amplitude  and  residual  frequency  modulation 
effects  of  the  oscillator.  Several  reasons  for  the 
residual  FM  and  the  steps  we  took  to  reduce  its  effect 
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are  presented.  Other  aspects  of  our  FM-CW  polarimetric 
radar  are  also  presented.  A  comparison  of  the 
performance  of  this  radar  and  a  step  frequency  radar  is 
given,  as  is  a  technique  for  determining  copolar  phase. 

Chapter  5  includes  a  discussion  of  analysis  of  ice 
surface  truth  data.  Included  are  computation  of 
autocorrelations  for  surface  and  volume  parameters. 
Simulations  are  done  to  illustrate  uncertainty  in  the 
autocorrelation  data.  Consideration  is  also  given  to 
statistics  of  the  inhomogeneities  in  the  desalinated 
ice . 

A  comparison  of  results  with  theory  is  the  subject 
of  Chapter  6.  Measured  surface  parameters  are  used  in 
the  small  perturbation  and  Kirchoff  (scalar 
approximation)  surface  scatter  models  to  generate 
confidence  intervals  with  which  to  compare  saline  ice 
backscatter  data.  The  results  and  a  simple  model  are 
used  to  show  that  the  volume  scatter  contribution  is 
important  for  very  smooth  saline  ice  surfaces.  The 
return  from  desalinated  ice  is  modelled  as  a  sum  of 
surface  and  volume  scatter  contributions.  Surface 
scatter  is  modelled  with  the  above-mentioned  models;  the 
volume  scatter  contribution  is  modelled  with  a  dense 
medium  theory.  Snow  cover  for  saline  ice  is  modelled  as 
in  Ulaby  et  al .  (19821. 

General  conclusions  and  recommendations  for  future 
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research  are  given  in  Chapter  7. 
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2.0  Background 

2.1  The  dielectric  properties  of  sea  ice  at  microwave 
frequencies 

An  extensive  literature  search  was  conducted  to 
survey  the  current  knowledge  of  the  electrical 
properties  of  sea  ice.  Because  so  much  work,  both 
theoretical  and  experimental,  has  been  done,  only  an 
overview  is  presented  here.  For  a  thorough  discussion 
of  the  subject  see  Vant  [1976]  and  Arcone  et  al  [1986b]. 

Because  the  electrical  properties  of  sea  ice  depend 
strongly  on  morphology,  we  first  give  a  description  of 
the  formation  and  morphology  of  sea  ice. 

2.1.1  Formation  and  morphology  of  sea  ice 

Sea  water  consists  primarily  of  H20,  Na+  and  Cl” 
with  smaller  percentages  of  SO*  ,  Mg2*,  Ca2+"  and  K+ .  The 
maximum  salinity  (in  grams  of  dissolved  solids/kg  of 
solution)  of  ocean  water  is  about  35o/oo  but  may  be 
lower,  depending  upon  dilution  with  fresh  water.  Figure 

2.1  indicates  the  change  in  the  sea  water  make-up  as  ice 
forms  to  equilibrium  at  given  temperatures.  The 
solution  consists  of  pure  ice  and  brine  down  to 
approximately  -24°C;  at  colder  temperatures  NaCl  (the 
principal  salt)  precipitates. 

If  sea  water  is  cooled  slowly  enough,  complete 
separation  of  ice  and  brine  can  occur.  This  rarely 
happens  in  nature,  especially  in  the  polar  regions. 
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Temperature  (C) 

Figure  2.1:  Phase  relations  for  ’’standard”  sea  ice. 

Circles  on  the  brine-salt  line  indicates 
temperatures  at  which  solid  salts 
precipitate  (Assur,  1958). 


where  ice  formation  occurs  rapidly  enough  to  allow 
considerable  amounts  of  brine  to  be  trapped  within  the 
ice  structure.  The  amount  trapped  depends  largely  on 
the  rate  of  ice  growth  CNakawo  and  Sinha,  1981 J;  the 
primary  entrapment  occurs  at  interplatelet  boundaries 
within  individual  crystals  composing  the  ice  structure. 

The  simplest  ice  structures  consist  of  two  layers  as 
shown  in  Fig.  2.2.  The  consistency  and  depth  of  the 
upper  layer  is  dictated  primarily  by  the  cooling  rate 
and  turbulence  in  the  water.  The  lower  layer 
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preferentially  grows  as  long  columnar  crystals,  each 
•  crystal  consisting  of  similarly  oriented  platelets 

(Martin,  1979]. 

When  the  sea  surface  is  calm  and  rapidly  cooled, 

^  the  water  crystallizes  in  small,  needle-like  shapes. 

These  rapidly  change  to  circular  disks  which  reach  a 
diameter  of  no  more  than  2  to  3  mm  [Weeks  and  Assur, 
1969].  Under  these  conditions  the  upper  layer  continues 
f  to  a  depth  of  1  to  2  cm. 

When  a  cold  wind  blows  over  sea  water  that  has  a 
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temperature  near  freezing,  the  water  is  stirred  and  ice 
crystals  of  1  to  10  Mm  in  thickness  and  1  to  2  mm  in 
diameter  form  in  the  water  column.  These  crystals  float 
to  the  surface  and  form  a  porous  mass  with  a  random 
crystal  structure  [Martin,  1979].  This  type  of  ice, 
referred  to  as  frazil,  can  form  to  thicknesses  of  up  to 


1  m . 


Ice  similar  in  structure  to  frazil  forms  when  snow 


falls  into  sea  water  that  is  supercooled.  Further 
cooling  causes  the  ice  crystals  formed  from  the  snow  to 
congeal  into  slush  ice.  This  type  of  ice  can  also  form 
layers  up  to  about  1  m  thick. 

The  locations  of  the  brine  inclusions  within  the 
ice  structure  are  illustrated  in  Fig.  2.2.  Typical 
brine-  pocket  dimensions  in  the  columnar  ice  are  on  the 
order  of  0.05  mm  [Anderson,  1960;  Pounder,  1965]  in 
diameter  and  2  to  3  cm  in  length  [Pounder,  1965; 

Addison,  1969].  Deep  within  the  ice  the  brine  pockets 
are  elongated  vertically. 

The  brine  content  of  sea  ice  is  not  static, 
however.  During  the  first  few  days  after  the  ice  forms, 
the  brine  content  of  the  ice  can  decrease  dramatically; 
this  decrease  is  a  f i rs t -generat i on  desalination 
process.  Ice  that  has  only  undergone  this  type  of 
desalination  is  referred  to  as  first-year  ice. 

Increased  temperatures  and  melting,  leading  to 
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formations  of  melt  ponds  during  the  spring  and  summer 

months,  cause  additional  reduction  in  the  ice  brine 
content  of  the  ice;  this  decrease  is  second-generation 
desalination.  Ice  that  has  undergone  at  least  one 
season  of  this  type  of  desalination  is  generally 
referred  to  as  multiyear  ice. 

First-generation  desalination  in  winter  ice  is 
illustrated  in  Fig.  2.3  -  notice  the  magnitude  of  the 
salinity  decrease  that  occurred  during  the  week  after 
the  ice  formed.  Four  mechanisms  to  explain  salinity 
decrease  i *•  ea  ice  have  been  presented  by  Untersteiner 
[1968]  1)  brine  pocket  diffusion,  in  which  a  salinity 

gradient  causes  continuous  preferred  melting  of  the 
brine  pocket  boundary;  (2)  gravity  drainage;  (3) 
flushing;  and  (4)  expulsion,  caused  by  pressure  build-up 
within  a  brine  pocket  as  cooling  occurs.  The  brine 
pocket  diffusion  hypothesis  is  not  considered 
significant  due  to  the  rapid  rate  at  which  desalination 
occurs.  The  gravity  drainage  explanation  is  believed  to 
be  a  significant  mechanism  over  the  long  term.  The 
flushing  mechanism,  due  primarily  to  snow  melt  at  the 
ice-snow  interface,  is  important  during  the  spring  and 
summer  but  is  not  effective  at  the  temperatures 
encountered  during  the  growth  of  winter  sea  ice  sheets. 
Expulsion  has  been  proposed  as  the  dominant  first- 
generation  desalination  mechanism  [Untersteiner,  1968]. 
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Figure  2.3:  Variation  of  salinity  with  time  for  2.5  cm 
sections  of  ice  at  a  depth  between  40  and 
42.5  cm.  Broken  line  indicates  "stable" 
salinity  or  average  of  all  the  salinities 
except  the  intial  high  value  at  this  level 
(Nakawo  and  Sinha,  1981). 


Expulsion  can  be  explained  as  follows:  further  cooling 
after  initial  brine  pocket  formation  causes  water  to 
freeze  out  of  the  trapped  brine  leaving  a  more 
concentrated  brine  and  a  smaller  pocket;  the  increased 
pressure  due  to  volume  reduction  forces  brine  out  of  the 
pocket . 

A  typical  salinity  profile  for  first-year  ice  is 
shown  in  Fig.  2.4  -  notice  that  the  curve  is  C-shaped. 
Brine  expulsion  from  very  thin  ice  causes  a  thin  (few 
tenths  of  a  mm)  high  salinity  (reported  to  be  as  high  as 
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Figure  2.4:  Salinity  and  temperature  profiles  of  a 
typical  frazil  ice  core  (Martin,  1979). 


95o/oo  [Martin,  1979])  layer  to  form  on  the  ice 
surface.  The  salinity  of  the  ice  within  a  few  cm  of  the 
ice-water  interface  is  influenced  largely  by  the  brine 
content  of  the  underlying  water.  The  salinity  of  the 
intermediate  structure  is  primarily  dependent  on  the 
growth  history  and,  apparently,  also  on  the  ice 
structure.  Nakawo  and  Sinha  [1981]  indicate  that  fine¬ 
grained  ice  retains  more  brine  than  columnar -gra ined  ice 
for  a  given  growth  history. 

First-year  ice  can  reach  a  thickness  of  up  to  about 
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0  2m  [Parashar,  19741;  before  it  attains  a  greater 

thickness,  the  ensuing  spring  and  summer  bring  about 
bulk  structural  and  salinity  changes.  Spring 
temperatures  cause  melting  at  the  ice  surface,  both  of 
the  snow  cover  and  of  the  ice  itself,  resulting  in  an 
altered  terrain.  The  warmer  temperatures  and  free 
surface  water  cause  the  brine  drainage  channels  to  grow 
®  and  expand  (See  Fig.  2.5).  These  are  present,  but 


poorly  developed,  in  cold  first-year  ice.  The  result  is 
that  the  meltwater  flushes  brine  from  the  ice.  This 
second-generation  desalination  results  in  a  salinity 
profile  typical  of  that  shown  in  Fig.  2.6.  Also  this 
flushing  of  brine  pockets,  and  entrapped  air  in  the  snow 
layer,  cause  air  bubbles  to  replace  brine  down  to  the 
freeboard  level.  The  resultant  ice,  ice  which  may 
persist  for  many  years,  is  referred  to  as  multiyear  ice, 
and  can  reach  thicknesses  of  several  meters. 

2.1.2  Electrical  properties  of  sea  ice 

Dielectric  measurements  have  been  made  at  a  variety 
of  frequencies,  salinities,  and  temperatures.  Vant  et 
al.,  [1978]  tabulated  the  conditions  under  which  many  of 
these  measurements  were  made  and  commented  on  the 
techniques  used  (See  Table  1.1).  Since  1978,  other 
experimenters  have  contributed:  Troy  [1981]  conducted 
dielectric  measurements  at  90  GHz  and  below;  Timco 
[1979]  studied  the  in-situ  resistivity  of  sea  ice;  and 
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a  Ice  dielectric  constant 
prior  to  1978  (Vant  et  al. 


Arcone  [1986a]  examined  the  dielectric  properties  of 

saline  ice  at  microwave  frequencies. 

In  the  frequency  range  of  1  -  18  GHz,  the  real  and 
imaginary  components  of  the  dielectric  constant 
generally  decrease  with  frequency.  This  is  primarily 
because  the  large  dielectric  constant  of  brine  decreases 
with  frequency  over  this  range  (See  Fig.  2.7). 
Hallikainen  [1982]  derived  a  set  of  multipliers  to 
describe  the  change  of  the  real  and  imaginary  components 
of  the  dielectric  constant  as  a  function  of  frequency: 
e'(18  GHz)=0.9  8«'(5  GHz )  = 

0 . 9  3  c'(  1  GHz  )  =0 . 8  6  c'(  0 . 5  GHz)  (2.1) 

1 8  GHz)=0.92t"(5  GHz  )  = 

0.61«"(1  GHz  )-0.46*'t0.5  GHZ)  (2.2) 

Vant  [1976]  determined  the  dielectric  constants  for 
many  cores  taken  during  the  Arctic  Ice  Dynamics 
Experiment  (AIDJEX)  and  found  a  decrease  in  both  real 
and  imaginary  components  as  a  function  of  temperature; 
his  measurements  were  at  1  -  4  GHz  and  salinities  of  4 
to  9o/oo.  The  real  component  decreased  from  about  3.7 
to  3.2,  and  the  imaginary  component  decreased  from  about 
0.4  to  0.03  over  the  range  of  -5°  to  -40°C.  These 
results  are  not  surprising  considering  that  the  amount 
of  liquid  in  the  brine  in  the  ice  varies  directly  with 
temperature . 

The  dielectric  constant  of  sea  ice  is  not 
isotropic.  Vant  et  al.  provide  many  curves  to 
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Figure  2.7:  Real  and  imaginary  parts  of  the  relative 

dielectric  constant  of  brine  (Vant  et  al., 
1978  )  . 


illustrate  the  effect  of  brine  pocket  orientation  on  the 
dielectric  constant,  as  predicted  by  theory  [Vant, 1976; 
and  Vant  et  al . ,  1978  ]  . 

The  dielectric  constant  of  sea  ice  has  been 
modeled,  both  empirically  and  theoretically.  Vant  et 
al.  attempted  to  include  the  available  knowledge  on  the 
subject  into  a  comprehensive  semi -empir i ca 1  model 
[Vant, 1976;  and  Vant  et  al.,  1978).  The  empirical 
components  of  the  model  relate  the  relaxation  time. 


static  dielectric  constant,  and  conductivity  of  the 


brine  to  temperature  and  frequency.  They  also  relate 

brine  volume  to  salinity,  temperature,  and  density.  The 
theoretical  components  consist  of  the  Debye  equations 
for  brine  and  ice,  to  determine  the  corresponding 
dielectric  constants,  and  a  mixture  equation.  The 
latter  comes  from  a  solution  to  Laplace's  equation  for 
an  arbitrary  spheroid  (quasi-static  approximation). 

This  comprehensive  model  performs  veil  for  FY  ice 
but  has  not  performed  well  to  date  on  MY  ice  due  to  the 
variable  density  of  MY  ice  and  also  because  the  shapes 
of  brine  inclusions  in  MY  ice  are  not  well  understood 
[Vant,  1976;  and  Vant  et  al .  ,  1978). 

This  model  works  veil  in  the  range  400  MHz  -  7.5 
GHz  but  is  lengthy  and  requires  a  detailed  knowledge  of 
the  ice  structure.  For  many  purposes  a  simple  empirical 
model  is  adequate.  The  form  of  the  model  Vant  proposed 
is 

t-  a  +  bvbr 

€=  c  +  dVbr  (2.3) 

where  a,b,c,  and  d  are  constants,  and  Vbr  is  the  brine 
volume.  Vbr,  can  be  determined  using  empirical 
equations  developed  by  Frankenstein  and  Garner  (1967): 
Vkf.  =  S  (  (  -  49 . 18  5/T )  +  0.532)  parts  per  thousand,  (2.4) 
where  S  and  T  are  measured  salinity  and  temperature. 

The  constants  a,b,c,  and  d  as  determined  by  Vant  are: 
a  =  3.05,  b  =  .0072,  c  =  .024,  d  =  .0033  at  4  GHz 


a  =  2.46,  b  =  .0224,  c  =  .006,  d  =  .01  at  10  GHz.  The 
dielectric  constants  at  other  frequencies  can  be 
estimated  from  equations  (2.1)  and  (2.2)  above. 


The  empirical  equations  are  used  throughout  this 


paper 


2.2  Theories  used  in  remote  sensing  of  sea  ice 


Many  surface-  and  volume-scattering  theories  have 
been  applied  to  the  problem  of  understanding  and 
interpreting  microwave  backscatter  from  sea  ice.  The 
most  commonly  used  surface  scattering  theories  are  the 
Kirchhoff  (scalar  approximation)  and  small  perturbation 
methods.  Both  analytical  wave  and  radiative-transfer 
methods  have  been  used  to  model  the  volume-scattering 
component . 

2.2.1  Review  of  surface  scattering  theories 


2 . 2 . 1 . 1  Kirchhoff  surface  scattering  theor 


The  Kirchhoff  scattering  theory  is  a  physical 
optics  approach;  i.e.,  it  is  based  on  vector  second 
Green's  theorem  which  states  that  the  scattered  field  at 
any  point  within  a  source-free  region  bounded  by  a 
closed  surface  can  be  expressed  in  terms  of  the 
tangential  fields  on  the  surface  (Ulaby  et  al.,  1982]. 
The  development  begins  with  (Pg.  926  of  [Ulaby  et  al.. 


1982  ]  )  : 


E  =Kf?jXf((  nxE-  ^3ns  x  (  nxH  )  ]  )exp(jksr-ns  )dS 


(2.5) 
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where  K  =  - jk, exp( - jk5 R0 )/4TTR0 

=  unit  vector  in  scattered  direction 
n=unit  vector  normal  to  the  interface  inside  the 
medium  in  which  scattering  is  being  considered 
*775  =  intrinsic  impedance  of  the  medium  in  which  Es 
is  evaluated 

k5 =wave  number  of  the  medium  in  which  E1  is 
evaluated  (c.g.j  ih  medium  1) 

Rc=range  from  the  center  of  the  illuminated  area 
to  the  point  of  observation 
E,H  =  total  electric  and  magnetic  fields  on  the 
interface 

Assume  a  plane  wave  of  the  form  E1 =aE0 exp ( - jkt n* •  r ) 
incident  on  the  surface  and  a  coordinate  system  which 
includes  local  surface  slopes;  e.g.,  t  =  nj  xn/  I  ( nj  xn )  I , 
d=nixt,  and  nj=txd. 

The  theory  is  developed  assuming  that  the  fields  at 
a  point  on  the  surface  can  be  expressed  as  the  sum  of 
incident  and  reflected  fields,  where  reflected  fields 
are  related  to  the  incident  fields  via  the  Fresnel 
reflection  coefficient.  This  implies  that  the  surface 
radius  of  curvature  is  sufficiently  large  (locally)  that 
the  Fresnel  reflection  coefficient  is  applicable. 

More  approximations  are  then  made  assuming  other 
physical  characteristics  of  the  surface.  For  surfaces 
that  have  small  slopes  a  scalar  approximat ion  is  used. 
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In  the  scalar  approximation  the  slope  terms  in  t 
and  d  are  ignored;  and  the  reflection  coefficient,  which 
is  dependent  on  surface  slope,  is  expanded  in  a  Taylor 
series  with  respect  to  the  slope  terms.  E  can  then  be 
expressed  in  the  form 

E  *=/  Up<,  exp(  jkx  ( ns  -?$£  )  t  (  r-r  '  )  )dS  '  (2.6) 

where  Up*,  is  of  the  form 

Up(,  =a0  +at  Z„+aj  Zy  (Pg.  937  of  [Ulaby  et  al.,  1982  ]). 

An  ensemble  average  of 

Es«  E*  results  in  a  series  of  terms,  the  first  of  which  is 


‘"-‘M 


LLL-y  L-*" 

[exp(jqx(x'-x")+jqy(y'-y")  ) 


-L-t-L-t 


x  <exp(jqz(z(x' ,y' )-z(x",y" ) )>]  dx ' dy ' dx"dy" 


ZL2L 


=  la_  f  ft  (  2L-  lu  I  )  (  2L-  |v  |  ) exp (  jqxu  +  jqy v ) 


JJ 


-ZI-2L 


x  exp(-q|  «xz(  l-/»)  )dudv  (2.7) 

where  qx,qy,qt  are  propagation  constants,  2L  x  2L  is  the 
illuminated  area,  and  p  is  the  correlation  function. 

Note  that  the  exp(-q|  az(l-^>))  term  was  obtained  from 
<exp( jqx( z ( x ' , y ' ) -z ( x" , y" ) ) >  assuming  that  z(x',y')  and 
z(x",y")  are  jointly  gaussian  with  variance  w2  and 
correlation  coefficient  p.  The  remaining  terms  of  Up<? 
are  developed  similarly;  these  are  functions  of  exp(- 
q2  <jx(1-^))  and  derivatives  of  p. 
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So  that  an  analytical  form  can  be  obtained,  the 

term  exp(-q£  »a(l-/»))  is  expressed  as  series  in  /*  su*_h 
that  the  first  intensity  term  is 


I~la0llA0exp(-q*  ax  )  £  (qi«r1)h/n! 

n  =  0 

2L2L 


•  =4TTRa  <  IE  ^  1%/Ao  lEo|X. 


x  JJZ*  expi  Dq*u+D  qyv  )auav  (2.8) 

-1L-2L 

where  A0=(2L)1  is  the  illuminated  area  (See  Ulaby  et 
al.,  (  1982  ]  for  development  of  other  terms).  After  a 
functional  form  is  assumed  for  p,  o-p^  is  determined  as 

(2.9) 

2. 2. 1.2  Small  Perturbation  method 

This  theory  begins  with  the  wave  equation  (Page  949 
of  (Ulaby  et  al.,  1982]),  i.e., 

VJE  +  k1E=0.  (2.10) 

Using  rectangular  coordinates,  each  component  is 
expressed  as  a  function  of  wavenumbers.  For  example, 


Ef)<  =(x/2TT)|'JExexp  ( -  jk.„x-jkyy)dxdy 


(2.11) 


(Ex  can  be  retrieved  from  E  f%  via  an  inverse  Fourier 
Transform)  Using  the  wave  equation  can  be  expressed 


as 


(-k*  -ky  +da/dz*  )Ejx +k2  Efx  =0, 
which  has  a  solution  of  the  form 
Aexp( - jk^z ) +A ' exp ( - jk*  z ) . 

The  total  fields  are  assumed  to  be  of  the  form 


(2.12) 
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Ex=(l/2TT)|‘j,Ux  exptjk^x  +  jkyy-jkfcZjdkxdky 

Ey=0-/2TTj|*JUy  exp( jk*  x  + jky  y- jkx z )dk*dky 

+exp( - jk*  sinS) (exp( jkxcos0)+Rexp( - jkzcos6) ) ) 
Ei=0l/2TT)|,JLJ2  exp(  jkxx+ jkyy- jkxz  )dkydky 


(2.13) 


where  Ux/Uy,Uz  represent  unknown  amplitude  with  respect 
to  k„  and  ky  and  the  additional  term  in  Ey  accounts  for 
the  source.  The  fields  in  medium  2  are  similarly 
expressed,  but  with  unknown  amplitudes  D*,Dy,  and  Dz 
(Notice  that  propagation  is  in  the  x,z  plane  such  that 
the  incident  wave  is  horizontally  polarized). 

From  the  boundary  conditions,  rix£E=0  and  nxdH=0, 
and  the  divergence  relations,  t7-E=0,  six  equations  are 
obtained.  AEx,£Ey,  and  £ Ez  can  be  determined  by 
expressing  UX,UV,UX,  D„,Dy,  and  Dx  in  terms  of  a 
perturbation  series  and  expanding  exp(-jkzz),  since  kxz 
is  small;  for  example, 

Ex-  (1/  UX1  +U*x  +  •  •  •  i  t 1- jkzz- .  .  .  ]exp(jkxx  +  jkyy)dkxdky 

Ex  =  Cl  /  2  TT^j* J*  C  DX1  +  DyZ  +  .  .  .  1  1 1-  jk*z-  .  .  .  lexp(jkxx+jkyy)dkxdky 


so  that 


^X-(1/2TT)|,,J*[UX1  -D>t  ]exp(jkxx+jkyy)dk*dky=0 
( to  first-order ) . 

First-order  terms  for  ^E  y  and  £E-z.are  obtained 


(2.14) 
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similarly.  AH*,£Hy/  and  &HZ  are  obtained  from 

derivatives  of  £E„,*^Ey,  and  AEZ  .  After  solving  the 
system  of  six  equations  and  six  unknowns  the  first-order 
scattered  field  can  be  expressed  in  terms  of  D* ,  Dy ,  Dx, 
Ux ,  Uy ,  and  Uz  .  An  ensemble  average  is  then  taken.  The 
result  is  of  the  form 

<EpqE*p«,  >=  l2k«7iCose*p<,  |2JJW(kx +ksin6,ky  )dkxdky  (2.15) 

where  ff  (  ky  ,  k  y  )=(i/2TT)j*j,p(  u,  v  )exp(  -  jkx  u- jky  v)dudv 

and  where  p(u,v)  is  the  correlation  function  of  the 


The  resulting 

<rp<>=8  I  ka  ^cosecosSjOtp,  |2W(k)t  +  ksinG,ky).  (2.16) 

(See  page  961  of  [Ulaby  et  al.,  1982]) 

2. 2. 1.3  Validity  conditions 

The  Kirchhoff  model  (scalar  approximation)  is  valid 
for  kl>6,  11>2.76<tA,  and  (  -  tzp"  (  0  )  ]  0 . 5  <  0.25  (Ulaby  et 
al.,  1982], 

where  <r  =  rms  height 

1  =  correlation  length 
P  =  correlation  function 
A  =  electromagnetic  wavelength 
k  =  wave  number 

The  small  perturbation  model  is  valid  for  k a  <  0.3, 
and  (V?)*/!  <  0.3  (Ulaby  et  al.,  1982]. 

2.2.2  Volume  scattering  theories 


One  of  the  most  commonly  used  volume  scattering 
theories  is  the  radiative-transfer  theory.  It  consists 
of  a  differential  equation  that  is  based  on  propagation 
of  intensity  and  is  not  suited  for  conditions  where 
fields  are  correlated. 

More  recent  emphasis  has  been  placed  on  analytical 
wave  theories.  These  can  generally  be  divided  into  two 
categories:  (1)  scattering  from  fluctuating  media,  and 
(2)  scattering  from  layers  containing  discrete 
scatterers . 

Fluctuat ing-medium  theories  include  the  Born 
Approximation  and  the  Renormalization  method.  The 
former  is  applicable  to  volumes  containing  small 
permittivity  fluctuations;  i.e.,  I  ( r ) /C^l < <1  where 
€  (  r* )  +  €  f  (  r  )  (Page  1067  of  [Ulaby  et  al.,  1986].  If 

the  electric  field  outside  the  medium  is  E  and  that 
inside  the  medium  is  Elf  then  VxVxE-kaE=0  outside  the 
medium  and  ^xVxEj^-  inside  the  medium. 

The  scattered  field  can  be  determined  by  iterative 
solution  of 

E  =  ECoJ  +jG0l(r,f1  )*f  (r^-E^r*.  >d? 
and 

E  =E[0)  +JGli  (z/z1)f^(rx  (2.17) 

(See  page  320  of  [Tsang  et  al.,  1985]) 
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where  G01(r,r1)  and  g41  (r,r1  )  are  layered-media  dyadic 

Green's  functions  in  the  absence  of  permittivity 
fluctuations  [Tsang  et  al.,  1985].  The  first-order  Born 
solution  (Born  Approximation)  is 


sl}(r)  *JG01(r,r1)ef  (r^-E^  (r* 


e  <05 


(2.18) 


)dr 

From  this  the  ensemble  average  of  E  )  is 

obtained  as  a  function  of  the  correlation  of  .  After 
a  functional  form  is  assumed  for  the  correlation 
function,  the  scattered  intensity  is  determined  and  <r° 
is  obtained. 

The  Born  Approximation  does  not  include  loss  due  to 
scattering  but  only  that  due  to  absorption.  When 
scattering  loss  must  be  considered  the  Renormalization 
method  is  used.  This  is  similar  to  the  Born 
Approximation  with  the  exception  that  the  attenuation 
due  to  scattering  loss  is  included  in  the  wave  equation; 
i.e.,  an  effective  propagation  constant  is  determined. 
The  equation  is  of  the  form: 

.  Z : 


VxVxEt  =  k  E1  =  -JC(r,r1  )G1±  (r,r1_)*E1(r1)dr 


(2.19) 


(See  page  1073  of  [Ulaby  et  al.,  1986)) 

where  C(r,r2)  is  the  correlation  of  the  dielectric 

constant  (the  integral  term  is  Dyson's  equation). 

Effective  medium  approximations  are  also  commonly 
used  for  layers  of  discrete  scatterers.  A  mixing 
formula  is  used  to  obtain  the  effective  dielectric 
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constant.  To  account  for  scattering  loss  ia  nation 

power  of  the  particles  is  considered.  Scattering  loss 
is  dependent  on  the  pa ir-d istr ibut i on  function,  just  as 
it  was  dependent  on  correlation  length  for  the 
Renormalization  method.  With  these  considerations  an 
effective  propagation  constant  for  spherical  particles 
can  be  found  to  be  [Tsang  et  al.,  19851: 
K2=k2(l+2fy)/(l-fy)+j  2  f  kraJ ly/(l-fy)  \z 
x  {  1  -8  f  +  2  4  f  H0  )  (2.20) 

where  f=n0v0  is  the  fractional  volume  occupied  by  the 
scatterers 

y=( €s-*)/( €s+20 ,  where  es  =  permittivity  of 

the  scatterers 
£=  permittivity  of 
the  background 

a=particle  radius 

H0  includes  pair  distribution  effects 
After  the  development  is  complete  <r°  of  the  form 
*■  v v  "%°„  =  I  (K-kiz_  )ki2  /n0 \Z  (n0/2nim(Kx)  ) 

(g 

x  { l  +  n0J*[  g  (  r  )  - 1  ]  dr  }  (2.21) 

-  co 

(See  page  438  of  [Tsang  et  al.,  19851) 
is  obtained  where 
k^2  =  kcos0£ 

K^=sqrt(Kz-k2sina©£  ) 
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g ( r ) =tvo-part icle  distribution  function 


More  recently  theories  have  been  developed  to 
account  for  multiple  scattering.  These  are  not 
considered  in  this  paper. 
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3.0  Measurements  conducted:  microwave  backscatter  and 

ice  characterization 

We  performed  backscatter  measurements  of 
controlled  and  we  11 -character ized  ice  surfaces  at  the 
Cold  Regions  Research  and  Engineering  Laboratory 
( CRREL ) ,  in  Hanover,  New  Hampshire.  CRREL  is  well 
suited  for  such  experiments  both  because  of  its 
excellent  facilities  (esp.  for  characterizing  ice)  and 
because  of  the  glaciological  expertise  available  on¬ 
site.  In  addition  the  local  climate  usually  favors  ice 
growth . 

Simulated  sea  ice  was  grown  at  CRREL  in  a  12.2  x 
5.2m,  1.5-m  deep,  outdoor  pond  that  can  be  covered  by  a 
large  moveable  tent  during  periods  of  inclement 
weather.  Straddling  the  pond  is  a  gantry  on  which  the 
radar  antennas  were  mounted,  permitting  collection  of 
data  over  the  entire  surface  at  various  incidence 
angles.  Adjacent  to  the  pond  is  a  tent  for  operating 
the  radars.  The  setup  for  the  1987  experiment  is  shown 
in  Fig.  3.1.  Figures  3.2  and  3.3  are  photographs  of 
setups  used  during  the  1988  and  1989  experiments, 
respectively. 

Table  3.1  provides  a  list  of  backscatter 
measurements  conducted  during  the  winters  of  1986-1987, 
1987-1988,  and  1988-1989.  During  all  three  campaigns 
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Figure  3.2:  Setup  for  CRREL'88  experiment.  Tent  (T) 

has  been  rolled  off  to  expose  pond  (P)  for 
measurements  from  gantry  (G). 


Figure  3.3:  Setup  for  CRREL'89  experiment.  University 
of  Kansas  and  Environmental  Research 
Institute  of  Michigan  radars  mounted  on 
the  gantry  (G),  University  of  Washington 
(UW)  radiometers,  and  University  of 
Massachusetts  (UM)  radar  and  radiometer 
during  measurements  on  pond  (P).  Weather 
station  (W)  is  in  foreground. 
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I&3I. 


Sals 


Scene 

Observed 


Comments 


• 

1987 

Jan.  17 

Jan.  18 

Smooth 
saline  ice 

Snow  covered  smooth 
saline  ice 

X-band 
angle  scan 

If 

Jan.  21 

Roughened  saline 

ice  " 

• 

1988 

Jan.  14,  15 

Smooth 
saline  ice 

C-band 
angle  scan 

Jan.  15 

Rough  saline  ice 

If 

• 

Feb.  26 

Desalinated  ice 

C-band  vs 
temperature 

Feb.  27,  28 

Desalinated  ice 

C-band 
angle  scan 

• 

1989 

Feb.  4,5 

Desalinated 
ice  blocks 

C-,Ku-band 
angle  scan 

Feb.  3 

Bare  and  snow  C-band 

covered  saline  ice  angle  scan 

• 

Table 

3.1  Scenes  observed  during  three 
experiments  at  CRREL 

years  of 

smooth  saline  ice  surfaces  (i.e.,  surfaces  with  rms 
height  <  0.05  cm)  were  observed;  thickness  effects  were 
stressed  during  1988  and  1989.  The  effect  of  a 
several-cm  snow  cover  was  studied  during  the  1987  and 
1989  campaigns.  The  effects  of  roughness  on 
0  backscatter  from  saline  ice  were  examined  in  1987  and 

1988.  Backscatter  from  desalinated  ice  was  the 
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emphasis  of  the  1988  and  1989  experiments. 


The  backscatter  measurements  of  1987  were  conducted 
at  X-band.  In  1988  we  concentrated  on  C-band  although  we 
made  some  X-band  measurements.  In  1989  very  thin  saline 
ice  was  studied  with  C-  and  Ku-  band,  thin  saline  ice  and 
snow  cover  effects  were  examined  with  C-band,  and  C-  and 
Ku-  band  were  used  to  study  backscatter  from  desalinated 
ice . 


The  ice  characterization  measurements  done  varied 
from  1987  to  1989  because  we  tried  to  provide  a  more 
detailed  and  more  accurate  quantitative  description  of 
the  ice  for  each  succeeding  experiment.  Discussions  of 
techniques  and  measurements  done  for  each  experiment 
are  given  in  the  following  sections.  For  all  of  the 
experiments  salinity  profiles  and  ice-structure 
photographs  were  supplied  by  CRREL .  Salinity  profiles 
were  determined  with  a  Beckmann  Salubridge  by  measuring 
the  conductance  of  the  meltwater  of  1-2  cm  thick 
sections  cut  from  an  ice  core.  Using  water  temperature 
as  an  input,  conversion  of  conductance  to  salinity  was 
performed  automatically  by  the  bridge.  Ice-structure 
photographs  were  obtained  of  0.5-1  mm  thick  ice 
sections,  which  were  prepared  using  a  Leitz  Wetzlar 
sledge  microtome.  The  ice  sections  were  photographed 
through  crossed  polarized  filters  so  as  to  enhance  the 
ice  crystalline  structure. 
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3.1  CRREL  1987 


3.1.1  Measurement  regime 

Experiments  began  on  January  17,  1987,  with 
backscatter  measurements  on  saline  ice  consisting  of  6 
cm  of  granular  slush  ice  overlaying  6  cm  of  congelation 
ice.  One  day  later,  after  6.5  cm  of  snow  had  fallen, 
further  measurements  were  made.  Three  days  later,  the 
third  set  of  measurements  involved  slightly  roughened 
ice.  Its  surface  differed  from  the  one  observed  earlier 
because  of  small-scale  roughness  introduced  by  the 
process  of  removing  all  snow  and  slush  with  a  rubber 
bladed  shovel. 

3.1.2  Experimental  procedure 

For  each  of  the  scenes,  backscatter  was  measured 

O 

at  incidence  angles  of  0,  10,  20,  30,  40,  and  50  for 
vertical  and  horizontal  polarizations.  The  height  of 
the  gantry  was  adjusted  to  maintain  a  slant  range  close 
to  3.5  m  for  all  incidence  angles.  Ten  statistically 
uncorrelated  observations  were  made  for  each  angle  by 
moving  the  gantry  over  the  ice  surface.  The  system 
response  was  determined  at  the  beginning  and  end  of 
each  of  the  three  measurement  sequences.  Echoes  from  a 
12-inch  Luneberg-lens  reflector  provided  an  absolute 
calibration. 

3.1.3  Measured  surface  physical  characteristics 

The  salinity  of  the  pond  water  was  maintained  at 
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24o/oo  throughout  the  experiment.  Air  temperature  was 
continuously  monitored  at  an  adjacent  weather  station. 
We  estimated  surface  height  and  measured  ice  thickness 
and  salinity  daily.  Ice  thickness,  salinity  profiles 
and  structural  Information  were  also  obtained  on  blocks 
of  ice  periodically  removed  from  the  pond. 

The  ice  sheet  was  grown  by  allowing  naturally 
precipitating  snow  to  accumulate  on  the  surface  of  the 
pond.  This  created  a  layer  of  slush  which  ultimately 
froze  to  form  saline  slush  ice.  Continued  freezing  of 
the  salty  water  directly  beneath  the  slush  ice  led  to 
the  formation  of  columnar-textured,  congelation  saline 
ice.  The  actual  progress  of  this  growth  is  illustrated 
in  the  series  of  thin-section  photographs  shown  in  Fig. 
3.4.  On  January  20,  1987,  the  ice  sheet  reached  a 
thickness  of  12  cm.  It  had  about  equal  parts  of  slush 
ice  and  congelation  ice.  The  slush  ice  was  composed  of 
two  layers;  coarse-grained  ice  with  platy  grains  up  to 
several  millimeters  in  diameter  overlaid  a  fine-grained 
component  with  grains  less  than  1  mm  in  diameter.  The 
photographs  show  that  some  modification  of  the  texture 
of  the  slush  ice  occurred  as  it  aged.  The  salinity 
characteristics  (Fig.  3.5)  and  crystalline  structure  of 
this  ice  sheet  closely  resemble  those  obtained  at  the 
earliest  stages  of  growth  of  Arctic  sea  ice  in  a 
slushed-over  lead  undergoing  quiet  freezing  (Gow,  1988). 
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A  Course  Grained 
B  Fine  Grained 
C  Columnar  Ice 


Sal  ini ty  X. 
8 


Figure  3.5:  Salinity  profiles  of  saline  ice  sheet 
grown  during  CRREL'87  experiment. 


Temperature  profiles  were  not  determined.  Using 
the  air  temperature  and  salinity  profiles,  we  estimated 
the  dielectric  constants  to  be  3.4-jO. 19,  3. 55-jO. 25, 
and  3. 35-jO. 17  for  January  17,  18,  and  21  respectively. 
Surface  roughness  was  estimated  by  observation;  we 
estimated  the  rms  height  to  be  0.1  cm  for  the  two  bare 
ice  sheets.  A  texture  was  noticeable  on  the  January  21 
surface . 


3.2  CRREL  1988 


3.2.1  Measurement  regime 

Experiments  began  on  January  11,1988,  with  C-  and 
X-band  backscatter  measurements  on  open  water.  They 
continued  as  the  ice  formed  to  a  thickness  of  about  1 
cm.  Measurements  were  repeated  for  ice  thicknesses  of 
3,  7,  and  10  cm.  Then  1  x  2  cm  ice  prisms  were  dumped 
on  the  sheet  to  simulate  rubble  field  roughness  in 
preparation  for  further  measurements. 

From  January  16  to  February  25  the  ice  was  allowed 
to  desalinate;  during  this  period  rain  and  meltwater 
contributed  to  the  flushing  of  salt  from  the  ice 
sheets.  Several  measurements  were  then  conducted  to 
study  the  effects  of  temperature.  These  were  followed 
by  incidence-angle  measurements. 

3.2.2  Experimental  procedure 

For  each  of  the  scenes,  backscatter  was  measured 
at  10,  20,  30,  40,  50,  and  60°  for  VV,  HH,  VH,  and  HV 
polarizations.  The  gantry  height  was  maintained  at 
approximately  3.3  m,  for  convenience.  For  each 
incidence  angle,  several  statistically  uncorrelated 
observations  were  made  by  moving  the  gantry  over  the 
ice  surface  -  the  number  of  uncorrelated  samples  varied 
depending  on  the  scene  we  were  observing.  We  measured 
the  system  response  after  every  2  or  3  samples. 
Backscatter  from  a  8"  Luneberg  lens  or  an  active  radar 


calibrator  was  measured  every  2  or  3  days. 
j.2.3  Measured  surface  physical  characteristics 

The  salinity  of  the  pond  water  was  maintained  at 
24o/oo  throughout  the  experiment,  as  during  1987.  Air 
temperature  was  continuously  monitored  at  a  nearby 
weather  station  and  readings  from  a  thermistor  string 
provided  ice  temperature  profiles.  We  measured 
surface-roughness  parameters  of  the  desalinated  ice  and 
obtained  saline-ice  roughness  parameters  from  other 
groups  participating  in  the  experiments.  We  also 
measured  ice  thickness  and  obtained  ice  salinity  and 
structural  information  from  CRREL . 

Ice  formation  began  as  a  result  of  spray  seeding 
of  the  surface  water.  The  internal  structure  of  the 
ice  sheet  after  about  6  days  growth  is  shown  in  Fig. 

3.6;  it  is  typical  of  that  of  first-year  ice  growing 
under  calm  conditions  and  composed  principaly  of 
columnar  crystals.  After  two  days  of  rather  intense 
desalination,  i.e.,  flushing  by  meltwater  and  rain,  the 
crystal  structure  changed  dramatically  and  more  bubbles 
appeared  (Fig.  3.7a).  Further  growth  and  desalination 
over  a  period  of  three  weeks  led  to  further  changes 
(Fig.  3.7b).  Two  and  a  half  weeks  later  the  ice  became 
very  porous  (Fig.  3.8). 

Platelet  structures  of  ice  crystals  as  observed  in 
horizontal  thln-sect ions  are  shown  in  Figs.  3.9  and  3.10 
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Figure  3.7:  (a)  saline  ice  sheet  of  Fig.  3.6  after  2 

days  of  intense  desalination  (b)  same 
sheet  after  3  weeks  of  growth  and 
desalination  (courtesy  of  CRREL ) . 
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Figure  3.8:  Porous  nature  of  ice  revealed  in  slice  of 
saline  ice  sheet  of  Fig.  3.6  photographed 
in  transmitted  light  after  5  and  a  half 
weeks  of  desalination  and  growth. 


Figure  3.9:  (a)  platelet  structure  of  crystals  at 

depth  of  6  cm  in  6.5-cm  thick  ice  sheet  of 
Fig.  3.6  after  3  days  of  growth  (b) 
platelet  structure  of  crystals  at  depth  of 
5  cm  in  14.5-cm  thick  ice  sheet  of  Fig. 

3.6  after  6  days  of  growth  (courtesy  of 
CRREL) . 
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Figure  3.10:  (a)  thermally  modified  platelet  structure 
in  ice  from  depth  of  8  cm  in  19-cm  thick 
ice  sheet  of  Fig.  3.6  after  11  days  of 
growth  and  desalination  (b)  thermally 
modified  platelet  structure  in  ice  from 
depth  of  8  cm  in  23-cm  thick  ice  sheet  ‘of 
Fig.  3.6  after  3  weeks  of  growth  and 
desalination  (courtesy  of  CRREL ) . 


The  platelet  boundaries  were  progressively  modified  as 
desalination  occurred.  Second-generation  desalination 
resulted  in  an  increase  in  air  bubbles.  The  desalination 
process  is  further  Illustrated  in  the  salinity  profiles 
of  Fig.  3.11. 

Figure  3.12  shows  temperature  profiles  obtained 
during  the  experiments.  Using  the  salinity  and 
temperature  profile  data  and  the  equations  of 
Frankenstein  and  Garner  [1967],  dielectric  profiles 
were  calculated  as  shown  in  Figs.  3.13.  Two-way 
absorption  loss  as-a-f unction-of -depth-into-the-ice  for 
each  ice  sheet,  obtained  from  the  estimated  dielectric 
constants,  is  shown  in  Fig.  3.14.  Table  3.2  presents 
surface  roughness  and  air  bubble  distribution  data. 

3.3  CRREL  1989 

3.3.1  Measurement  regime 

Experiments  began  on  January  9,  1989,  with  C-  and 
X-band  backscatter  measurements  on  open  water. 

Continuous  measurements  were  conducted  as  the  ice 
formed  to  a  thickness  of  2.5  cm.  Measurements  of  this 
smooth  ice  sheet  were  repeated  periodically  throughout 
the  remainder  of  the  experiment. 

The  major  focus  of  the  CRREL  '89  experiment  was  on 
microwave  measurements  of  desalinated  ice.  On  29 
December  1988  a  16-18  cm  thick  saline  ice  sheet  was  cut 
into  1  x  1  m  blocks  and  harvested.  These  blocks  were 
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real  epsilon 


(a) 


imag  epsi Ion 


(b) 


Figure  3.13:  Estimated  dielectric  constant  profiles  at 
C  band  on  dates  measurements  were  taken 
during  the  CRREL'88  experiment:  (a)  real 
component  of  epsilon  (b)  imaginary 
component  of  epsilon. 
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Figure  3.14:  Two-way  absorption  loss  at  C  band  as  a 

function  of  depth  into  the  ice,  computed 
from  data  of  Fig.  3.13. 


• 

surface _ 

rms 

he iaht 

s  .d . 
of  rms 
heiqht 

corr 

lenoth 

corr 

+  1 

s .  d  . 

length 
-1 
s  ,d . 

Smooth  saline 

0.029 

0 . 01 

1 . 77 

0 . 81 

2  .  49 

ice  ( Jan .  14 ) 
Roughened  saline 

0.45 

0 .14 

0.98 

0.56 

1.96 

• 

ice  (Jan.  15) 
Desalinated  ice 

0.18 

0.07 

2.5 

1 .88 

4 . 52 

(Feb.  26,27,28) 

volume 

mean  air 
bubble  radius 

air  bubble 

corr  lenoth  density 

Desalinated  ice  est.  unknown  unknown 

(Feb.  26,27,28)  67%  *0 . 5  mm 

32%  0.5-1  mm 
1%  2  mm 


Table  3.2  Surface  characteristics  for  CRREL  '88 
experiment  (all  dimension  in  cm  unless 
otherwise  indicated) 
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stored  on  pallets  and  allowed  to  desalinate.  These 
blocks  were  then  assembled,  as  shown  in  Fig.  3.15,  to 
allow  backscatter  measurements.  We  performed 
measurements  at  C-  and  Ku-  band  for  one  and  two  layers 
of  blocks. 

3.3.2  Experimental  procedure 

For  saline  ice,  backscatter  was  measured  at  10, 
20,  30,  40,  and  50°  for  VV,  HH,  VH,  and  HV 
polarizations.  For  the  desalinated  ice  block 
assemblage,  measurements  were  made  at  incidence  angles 
of  15,  25,  and,  for  Ku-band,  35°.  For  saline  ice,  5 
statistically  uncorrelated  observations  were  made  at 
each  incidence  angle.  Only  3  statistically 
uncorrelated  samples  were  possible  for  the  desalinated 
ice  because  we  performed  measurements  so  as  to  avoid 
the  effects  of  block  edges.  The  system  response  was 
determined  after  every  2  or  3  samples  and  backscatter 
from  a  8"  Luneberg  lens  provided  a  calibration 
reference . 

3.3.3  Measured  surface  physical  characteristics 

The  salinity  of  the  pond  water  was  maintained  at 
24o/oo  throughout  the  experiment,  as  in  1987  and  1988. 
Physical  characteristics  of  the  ice  were  measured  in 
nearly  the  same  manner  as  in  1988.  In  1989,  however, 
surface  roughness  was  determined  with  a  wire  gauge 
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Figure  3.15:  Desalinated  ice  blocks  as  studied  during 
the  CRREL’89  experiment. 


Figure  3.16:  Porous  structure  of  ice  as  seen  in 

horizontal  thick  section  from  top  surface 
of  a  desalinated  block  shown  in  Fig. 

3.15. 
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rather  than  from  photographs,  and  we  made  a  more 
thorough  analysis  of  the  air  bubble  structure  of 
desalinated  ice. 

Figure  3.16  illustrates  the  porous  structure  of 
desalinated  ice  as  observed  in  a  horizontal  section 
from  near  the  surface  of  a  block  and  Fig.  3.17  shows 
the  air  bubble  size  distribution  at  the  top,  middle  and 
bottom  of  the  block.  Other  structural  parameters  are 
listed  in  Table  3.3. 


surface 

rms 

he iaht 

s .  d . 
of  rms 
he iaht 

corr 

lenath 

corr 

+  1 

s  .  d  . 

length 

-1 

s  .  d . 

Bare  smooth  and 
snow  covered  saline 
ice  (Feb  3) 

0.024 

0.012 

1.2 

1.1 

1.46 

Desalinated  ice 
(Feb. 4  ,  Feb.  5) 

0.08 

0.15 

3.3 

2.2 

3.9 

dielectric 

air  constant 

volume _ bubble  radius  corr  length density 

Desalinated  ice  (see  Fig.  3.17)  1.5  -  2.0  mm  0.844 

(Feb.  4,  Feb.  5)  (bulk) 


mean 

Snow  Darticle  size 

wetness 

Feb.  3  . 1  x  . 2  cm 

0.05  g/cm3 

<5% 

sheet 

surface  rms  height  =  0.28  cm,  corr  length  =  1.8  cm 


Table  3.3  Surface  parameters  for  CRREL  '39  experiment 
(all  dimensions  in  cm  unless  otherwise 
indicated ) 
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Figure  3.18  indicates  the  temperature  profile  of 
the  saline  ice  sheet  at  the  time  of  the  backscatter 
measurements  and  Fig.  3.19  the  ice  sheet  salinity 
profiles  for  the  surfaces  observed.  We  obtained  the 
dielectric  constant  profiles  in  Fig.  3.20  using  the 
empirical  equations  given  by  Garner  and  Frankenstein 
[1967].  Two-way  absorption  loss  as-a-f unction-o£- 
depth- into-the-ice  for  each  ice  sheet,  obtained  from 
the  estimated  dielectric  constants,  is  shown  in  Fig. 


IgmperaturQ  (deg  C) 


Figure  3.18:  Temperature  profile  of  saline  ice  sheet 
grown  during  the  CRREL'89  experiment,  on 
date  measurements  were  taken. 


Figure  3.19:  Salinity  profiles  of  saline  ice  sheet  and 
desalinated  ice  blocks  on  date 
measurements  were  taken  during  the 
CRREL'89  experiment. 
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at  C  band) 

—  Feb  4  and  5 
'blocks  at  C 
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—  Feb  4  and  5 
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Feb  3.  1989 
(  saline  ice 
at  C  band) 
....  Feb  4  and  5 
(blocks  at  C 
band) 

—  Feb  4  and  5 
(blocks  at  Ku 
band) 


imag  epsilon 


(b) 


Figure  3.20:  Estimated  dielectric  constant  profiles  on 
dates  measurements  were  taken  during  the 
CRREL'89  experiment:  (a)  real  component 
of  epsilon  (b)  imaginary  component  of 
epsilon. 


67 


Figure  3.21:  Two-way  absorption  loss  as  a  function  of 
depth  into  the  ice,  computed  from  data  of 
Fig.  3.20. 
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4.0  The  Radar  Systems 


Over  the  course  of  three  years  of  experiments  at 
CRREL  three  different  radar  systems  were  used.  During 
the  1986-1987  experiment  the  X-band  SOURCESCAT  radar 
provided  backscatter  data;  during  the  1987-1988 
experiment  the  first  version  of  PAIR  (Polar imetric 
Agile  Instrumentation  Radar)  was  used;  and  during  the 
1988-1989  campaign  the  second  version  of  PAIR  and  a 
newly  developed  step  frequency  radar  (STEFRAD)  were 
employed . 

4.1  Design  considerations  and  operation 

The  constraints  on  the  designs  of  our  systems  were 
imposed  both  by  logistics  and  by  our  experiment 
objectives.  The  CRREL  facility  limited  measurement 
ranges  to  a  maximum  of  4  to  5  meters.  Our  experiment 
objectives  required  fine  resolution  both  in  slant  range 
and  in  azimuth.  The  1987-1988  and  1988-1989 
experiments  further  required  copolar  phase  measurement 
capability. 

4.1.1  The  SOURCESCAT  radar 

A  block  diagram  of  the  X-band  SOURCESCAT  radar  is 
shown  in  Fig.  4.1.  This  radar  is  a  simple  single¬ 
antenna  triangularly  modulated  FM-CW  device  ([Gogineni, 
1984]  and  (Zoughi,  1984]).  Fine  resolution  in  slant 
range  is  obtained  by  modulating  over  an  r.f.  bandwidth 
of  2  GHz,  giving  an  unweighted  range  resolution  of 
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Figure  4.1:  Block  diagram  of  X-band  SOURCESCAT  radar 
( Zoughi ,  1984). 


70 


Frequency 
r.f.  bandwidth 
FM  race 

range  resolution 
(weighted ) 
footprint 
(3.5  meter  range) 


10  GHz 
2  GHz 
100  Hz 
12  cm 

20  x  20  cm 


Table  4.1:  SOURCESCAT  radar  specifications. 


r res  =c/2B-7.5  cm.  Fine  resolution  in  azimuth  and  the 
plane  wave  condition  are  satisfied  at  ranges  within  the 
near  and  intermediate  fields  of  the  antenna  by  using  a 
parabolic  dish  antenna  with  feed  slightly  displaced 
from  the  focal  point.  Pertinent  system  specifications 
are  found  in  Table  4.1  (See  (Zoughi,  1984]  for  a  more 
detailed  discussion  of  this  radar). 

4.1.2  PAIR 

PAIR  operates  at  C  and  X  band  and  is  shown  in 
block  diagram  form  in  Fig.  4.2.  Notice  that  this  radar 
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X 


> 


a 


2 


Block  diagram  of  the  C-band  section  of 
PAIR  (the  X-band  section  is  nearly 
identical). 


operating  frequency 
r.f.  bandwidth 
modulation 
modulation  rate 
range  resolution 
(without  weighting) 
footprint 
(3.5  m  range) 
dynamic  range 
transmit  power 


6  GHz 
500  MHz 
triangular 
70  Hz 
30  cm 

20  x  20  cm 

>50  dB 
10  mW 


11  GHz 
500  MHz 
triangular 
70  Hz 
30  cm 

20  x  20  cm 

>50  dB 
10  mW 


Table  4.2:  PAIR  specifications. 


is  a  modified  and  enhanced  version  of  SOURCESCAT. 
Components  of  PAIR  not  found  in  SOURCESCAT  include: 

•  (1)  a  phase-locked-loop  for  frequency  linearization  and 
phase  stabilization;  (2)  I  and  Q  channels  to  provide 
amplitude  and  phase  information;  (3)  digital  data 

^  acquisition  system  to  provide  flexibility  in  signal 

processing;  and  (4)  mul t i po lar i za t i on  capability. 

Other  enhancements  include  a  focussed  offset-fed 
parabolic-dish  antenna  to  improve  side-lobe  and  cross- 

•  polarization  performance  (version  2  of  PAIR  only). 
Pertinent  system  specifications  are  found  in  Table  4.2. 
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4. 1.2.1  Discussion  of  circuitry  unique  to  PAIR 

A  phase-locked  loop  ( PLL )  is  used  to  improve  the 
performance  of  the  YIG  VCO .  It  does  this  by 
linearizing  the  f requency-ver sus -dr iving  current 
characteristic  of  the  oscillator  -  Caputi  [1973]  used 
this  technique  and  provided  an  analysis  of  this 
technique  for  Backward  Wave  Oscillators.  Oscillator 
nonlinearities  increase  the  amplitude  of  range 
sidelobes  and  the  change  in  such  nonlinearities  versus 
time  and  temperature  adds  uncertainty  to  phase 
measurements.  Figure  4.3  illustrates  the  frequency- 
versus-dr iving  current  ( f  vs .  i)  response  for  a  typical 
YIG  oscillator.  Maximum  At  in  Fig.  4.3  is  around  5-10 
MHz  for  the  oscillators  we  used.  Deviations  from 
linear  are  due  to  loading  of  the  YIG  phere  by  the 
active  element  and  to  interaction  between  resonance 
modes  within  the  sphere  [Clark  and  Swartz,  1971]. 

Figure  4.4  illustrates  the  effect  of  bias  level  on  the 
impedance  of  a  typical  active  element  (X-band 
oscillator  diode)  -  this  is  important  because  the  bias 
level  can  change  with  time  and  temperature.  This 
impedance  variation  is  sufficient  to  cause  output 
frequency  variations  of  2  to  3  MHz  per  volt  change  of 
bias  (Clark  and  Swartz,  1971].  Figure  4.5  shows  what 
happens  when  a  spurious  resonance  is  excited  within  the 
sphere  that  crosses  the  main  resonance  at  about  8.5  GHz 


FREQUENCY  -  GHz 


Figure  4.5:  Effect  of  a  spurious  YlG-sphere  resonance 
on  linearity  of  a  YIG-tuned  oscillator 
(Clark  and  Swartz,  1971). 


•  Figure  4.6:  Effect  of  temperature  on  linearity  of  a 

typcial  heated  YIG-tuned  oscillator  (Clark 
and  Swartz,  1971 ) . 
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[Clark  and  Swartz,  19711.  Figure  4.6  indicates  several 
MHz  variation  in  deviation  from  linear  over  the 
temperature  range  of  -55°C  to  100°C  even  when  the  YIG 
heater  is  used.  See  Clark  and  Swartz  (19711  for  a 
thorough  disscussion  of  factors  which  affect  YIG 
oscillator  linearity. 

Both  in-phase  (I)  and  quadrature  (Q)  mixer 
products  must  be  available  to  obtain  phase  information 
for  a  homodyne  receiver  such  as  ours.  A  commercially 
available  quadrature  mixer  with  phase  imbalance  of  less 
than  6°  and  amplitude  imbalance  of  less  than  0.3  dB  was 
selected  to  provide  I  and  Q  channel  outputs. 

Amplifiers  and  filters  for  the  two  channels  were 
carefully  designed  and  selected  so  that  balance  is 
maintained  within  one  or  two  tenths  of  a  dB  in 
amplitude  and  a  few  degrees  in  phase.  The  low-noise 
amplifiers  consist  of  two  stages.  The  first  stage  uses 
LM194  matched  transistor  pairs,  and  the  second  stage 
uses  AD625  instrumentation  amplifiers.  The  two  passive 
10-section  Bessel  filters,  a  2-kHz  low-pass  and  a  15- 
kHz  high-pass  for  each  channel,  are  commercial  units. 

Processed  I  and  Q  channel  outputs  are  digitized  to 
12  bits  at  a  40kHz  rate  by  an  Analog  Devices  RTI-860 
board.  A  Compaq  Portable  II  then  reads  the  data  from 
the  RTI-860  and  stores  them  on  floppy  disks  or  on 
Bernoulli  cartridges.  Programs  stored  on  the  Compaq 
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permit  computer  controlled  acquisition  of 
multipolar ized  data,  data  previewing  with  flexible 
windowing,  and  preliminary  data  analysis. 

Multipolarization  capability  was  added  because  of 
the  current  emphasis  on  polarimetric  measurements.  Use 
of  a  dual-polarization  feed  in  an  offset-fed  parabolic 
dish  conf igurat io i  provides  better  performance  than  a 
standard  center-fed  configuration:  higher  cross¬ 
polarization  sensitivity  and  lower  side-lobes.  A  dual¬ 
polarization  conical  feed  optimizes  the  antenna  cross¬ 
polarization  and  side-lobe  performance.  It  was  used 
Instead  of  a  diagonal  horn  because  of  the  convenience 
in  machining. 

4.1.3  STEFRAD 

To  obtain  better  performance  than  is  possible  with 
PAIR,  we  designed  and  built  a  step  frequency  rada’’  co 
operate  at  5.5  and  13.9  GHz.  Improved  stability  -s 
possible  with  this  type  of  radar  because  a  synthesized 
source  (i.e.,  a  source  phase-locked  to  a  crystal 
oscillator)  can  be  used.  Thus,  problems  associated 
with  linearity  and  phase  noise  can  be  greatly  reduced. 

By  incorporating  commercially  available  network 
analyzers  the  design  and  construction  of  such  radars 
becomes  relatively  easy.  A  high  degree  of  accuracy  in 
microwave  vector  measurements  and  dynamic  range 
exceeding  100  dB  are  additional  incentives  for  using 
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operating  frequency 

5.3 

GHz 

13.6 

GHz 

r.f.  bandwidth 

500 

MHz 

500 

MHz 

range  resolution 

30 

cm 

30 

cm 

(without  weighting) 
footprint 

1.2  x 

1.2  m 

1.2  x 

i .  2 

(3.5  m  range) 
dynamic  range 

>70 

dB 

>70 

dB 

transmit  power 

10 

raW 

10 

mW 

Table  4.3:  STEFRAD  specif ications . 


these  devices. 

A  block  diagram  for  STEFRAD  is  shown  in  Fig.  4.7. 
An  HP-8753A  network  analyzer  serves  as  the  core.  This 
unit  compares  the  target  return,  in  amplitude  and 
phase,  to  a  stable  reference  signal.  The  front-end 
consists  of  an  up-converter,  a  down-converter,  and 
appropriate  antennas  and  switching.  A  filter  in  each 
section  removes  unwanted  sidebands  and  an  amplifier  in 
each  section  compensates  for  mixing  and  cabling  losses. 
Horn  antennas  were  used  because  the  systems  have 
recently  been  built  and  are  just  now  moving  past  the 
initial  testing  stage.  Although  these  antennas  have 
wide  beamwidth  and  moderate  side-lobe  performance,  they 
perform  well  in  a  cross-polarized  configuration.  See 
Table  4.3  for  a  list  of  pertinent  system  parameters. 


Figure  4.7:  Block  diagram  of  the  Ku-band  section  of 
STEFRAD  (the  C-band  section  is  similar) 


4.2  Comparison  of  modes  of  operation 


4.2.1  SOURCESCAT 

The  SOURCESCAT  radar  is  triangularly  modulated. 

The  frequency  of  the  output  of  this  radar  is  thus 
proportional  to  target  range.  Consider  the  return  for 
a  point  target:  Let  the  transmitted  signal  be 
Acos(2TT£ct  +  2TTAft2  +  *)  .  (4.1) 

At  time  t  +  T  the  signal  becomes 

Acos  (  2TT£ c  ( t  +T)  +  2TTA£  ( t+T)2  +  *)  .  (4.2) 

Then  the  return  from  the  target  is 

Atcos[2n£t  (t  +  T)+2TTA£t(t+T)+*-2pczt-PsZs+**]  ,  (4.3) 

vhere 

At  =  amplitude  of  the  return  signal 

Af  =  time  slope  of  frequency  increase 
=  two-way  propagation  time 
£4  =  phase  shift  due  to  the  target 
*»c  =  2TT£c/c 

fc  =  lowest  r.f.  frequency  transmitted 
z^.  =  slant  range  to  target 

ps  =  effective  propagation  constant  within  system 
components  at  frequency  fc 
zs  =  effective  length  of  system  paths 
The  signal  at  the  radar  output  is  then 
(  AA*/2  )cos  (  (  2TTAfT)  t  +  2TTA£ T2+2/*c  z<  +  **y  zs  -  )  .  (4.4) 

If  Af =B/ ( Tm/2 ) , 
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where  B=r . f .  bandwidth 

Tm  =  modulation  period, 
and  Y=2R/c, 

where  R  is  the  range  to  the  target 
and  c  is  the  speed  of  light, 
then  the  frequency  of  the  radar  output  signal  is 
AfT=(B/( T„/2) ) (2R/C)  =  4BR/cT„.  (4.5) 

The  return  is  not  a  single  frequency,  however,  due 
to  the  change  in  phase  of  the  return  as  the  modulation 
slope  changes  sign.  The  actual  voltage  return 
spectrum,  as  determined  by  a  scalar  spectrum  analyzer, 
is  approximately 


F  ( f ) 


-lift 


,/2 


exp(  j  (  2TTA£T-tt»)t )  dt 


IJ?: 


(W2  exp(  j  (2TTAfT-w)t)  dt  I  )  x  ( f -nf  *,)  ) 

=  2(sin(  (  2TT£fT-ou)  (T„/4  )  )/(  2TTA£-W)  )x( *(  £-nf~)  ) , 

n=l, 2,3,...  (4.6) 
For  this  radar  we  obtained  data  using  an  HP3585A 
spectrum  analyzer  and  stored  the  data  on  mini  data- 
cassettes  via  an  HP-75  field  computer.  To  obtain  an 
accurate  estimate  of  the  target  return  envelope  we  set 
the  resolution  bandwidth  of  the  spectrum  analyzer  to 
include  at  least  one  spectral  line  (>fm)  and  set  the 
video  bandwidth  to  average  about  ten  up  and  ten  down 


sweeps  . 

Although  SOURCESCAT  is  a  simple,  reliable,  and 


easy-to-use  Instrument,  It  does  not  provide  for 


measurements  of  phase  and  flexibility  in  data  analysis. 
These  are  the  primary  reasons  that  PAIR  was  built. 

4.2.2  PAIR 

With  PAIR,  target  amplitude  and  phase  spectra  are 
obtained  by  performing  a  complex  FFT  on  the  average  of 
ten  time  samples  of  I  and  Q  channel  data  -  the 
averaging  is  done  to  reduce  system  effects.  The 
resulting  spectrum  for  a  point  target  is 
N-l 

F( k ) ( AAt/2  Jexpt  j  (  2T!Af  TnTs  +2TTAf  T*+2PcZ*  +  *szs-0*)  1 
n=0 

x  exp  (  -  j  (  2TTk/NTs  )  nTs  ) 

=  (AA4-/2)exp[  j  (2n^f  r2+2Ptz^  +  #s2I-0t  )  ] 
x  exp[  -  j  (  (  (N-l)/2)  (2TTAfr-(2nk/NTs  )  )TS) 1 
x  sinl  (NT,/2)  (2TTAfT-<2TTk/NTs  )  )  ] 

/  (sint  (Ts/2)  (2nAfT-(2TTk/NTs  )  )  I  (4.7) 

Target  amplitude  information  is  easily  obtained. 

Several  other  factors  must  be  determined,  however, 
before  the  target  phase,  0t,  can  be  extracted. 

Typically  we  are  most  interested  in  the  copolar 
phase  (0tvv-0tHH  Assuming  that  the  system  is  stable, 
measurement  of  a  known  target  such  as  a  sphere  allows 
us  to  find  the  target  copolar  phase.  For  a  sphere  we 
expect  0Vv=^hh-  Thus,  measuring  the  return  from  a 
sphere,  we  obtain 
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”^HHspK  -2TT^  Tvv  *"  Tmm  )+fs  (ZJw  _z;hm  ^ 

-(  (N-l  )/<..  (2n^f(TvvayJ  -Tn-sys  )  )  T*  (4.8) 

where  Ts  =  sampling  period 

Twsys  /  THHiys  =  system  delay  for  VV  and  HH 

channels 

Tw/  Thm  =  system  delay  plus  propagation  delay  to 
target  and  back  for  VV  and  HH 
polarizations 

2,Vv  /  zSHH  =  system  path  lengths  for  vv  and  HH 
channels 

For  a  stable  system  0s(z,vv-zSHM  )  and  2TTAf  ('■ftv,y5  -Th«*v»  ) 
will  be  nearly  the  same  when  the  target  measurement  is 
made.  For  short  ranges  and  values  of  Af  that  were  used 
the  2TTAf  ( Tv*  ~Tm*m  )  term  is  negligible.  Thus 

^  *  (4.9) 

For  this  type  of  scheme  to  produce  accurate  phase 
measurements  it  is  essential  that  the  system  be  very 
stable.  The  oscillator  must  be  stable  from  VV 
measurement  to  HH  measurement.  A  drift  of  5  MHz  in  fc 

Q 

causes  13  phase  difference  at  a  range  of  5  m. 

Stability  from  measurement  to  calibration  is  also 
important  but  not  as  critical  if  care  is  taken  in 
system  design.  A  drift  of  5  MHz  in  fc  in  this  instance 
only  introduces  an  error  of  about  2°  if  there  is  a  10 
cm  difference  between  the  HH  and  VV  channel  path 
lengths.  Ideally  one  should  phase-match  the  VV  and  HH 
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channels.  In  that  case  It  Is  not  necessary  to  use  a 

sphere  to  determine  the  copolar  phase. 

4.2.3  STEFRAD 

The  step  frequency  radar  operates  in  a  mode 
similar  to  that  of  PAIR.  At  short  ranges,  the 
resulting  equations  for  a  point  target  return  are 
identical.  Conceptually  the  only  difference  is  that  a 
step  frequency  radar  measures  return  amplitude  and 
phase  at  a  fixed  number  of  discrete  frequencies  whereas 
an  FM-CW  radar  measures  amplitude  and  phase  of  a 
frequency  shift. 

Iizuka  [1984]  explained  the  operating  mode  in 
terms  of  a  vector  addition  of  returns  from  several 
targets.  Consider  a  system  using  frequencies  of 
fcfnAf,  n  =  0, 1, 2,  . . . ,N-1 .  Measured  amplitude  and  phase 
at  fc+mAf  for  a  scene  consisting  of  four  targets  (See 
Fig.  4.8)  is 

Hm=soE0 / zj  exp(  j  (  2  (  2TT)  (  f t+mAf ) /c ) zD )  + 

st  E0/z*  exp  (  j  (  2  (  2TT)  (  f  c  +mAf  )  /c  )  zt)  + 
s2  E0/Zj  exp  (  j  (  2  (  2TT)  (  f  c  +mAf  )  /  c  )  z  )  + 
s3E0/z*  exp( j ( 2( 2n) ( f c+mAf )/c)z3 )  (4.10) 

where  s0,s1,s2,s3  =  voltage  scattering  coefficients 

of  targets  0,1,2,  and  3 

zo/Z1,z2,z3  =  slant  ranges  to  targets  0,1,2,  and  3 
E o  =  transmitted  signal 
For  the  case  of  N  targets 
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Figure  4.8:  Vector  addition  of  returns  from  four 
targets  (Iizuka,  1984). 


Hm=£  skE0/z£  exp  (  j  (  2  (  2TT)  (  f  c  +  n£f  ) /c  )  zk ,  («.ll) 

k  =  0 

Consider  the  case  of  N  targets  at  ranges  of 
z0,z»+Az,z0+2iiz,  .  .  . , z„+ (N-l ) Az  then 

H'n  =  (Hnexp(-j2(2TT)  (fc  +•  ndf)/c)zo)  =  [expC-j<2  (ZtiX-fci- 

N-l 

E( s„ Ee/z*  )exp[  j  (  2  (  2TT)  (  ft  +n£f  )  /c  )k  Az  ] 


x 


(4.12) 


If  we  let 


hK=(NskE9/z*  )exp(  j  (  2  (  2TT)  ( Ic/c  )k  Az  ]  (- 

then  H and  h ^  are  related  by  a  discrete  Fourier 
Transform.  That  is, 


(4.13) 


H’n=(l/N)  £  hkexp(j(2TT)(kn/N)]  and 
k=0 


hk  =  £  H'n  exp(-j(2TT)  (kn/N)}  (4.14) 

n=0 

where  N=c/2AfAz. 

Figures  4.9  -  4.11  are  used  to  further  illustrate 
how  a  step  frequency  radar  works,  assuming  8  stepped 
frequencies.  The  8  waveforms  in  Fig.  4.9  correspond  to 
phase  shift  encountered  by  the  r.f.  signal,  of 
frequency  fc+n£f,  in  propagating  to  the  target  and 
back,  relative  to  phase  shift  similarly  encountered  by 
the  r.f.  signal  of  frequency  f c ,  for  each  of  the  8 
stepped  frequencies  and  at  8  ranges  from  0  up  to  the 
maximum  unambiguous  range.  Note  that  one  period  of  the 
waveform  corresponds  to  a  2TT  phase  shift.  For  example, 
IT/4  phase  shift  occurs  at  £f  for  a  range  of  Ru/8,  17/ 2 
occurs  at  £>f  for  a  range  of  R^/4,  etc. 

The  vector  being  measured  by  the  network  analyzer 
is  of  the  form  Aexp(j$).  Consider  the  measurements  Ae , 
Atexp(  j  4TTAf  R  )  ,  Aaexp  (  j  8TTAf  R  )  ,  ...  ,  A?exp  (  j  28n&f  R  )  by 

an  8-frequency  step  frequency  radar  and  assume  that 
A0 =A^= . . . =A .  At  a  range  of  R=Ra/8,  the  measurements 
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Figure  4.9:  Relative  phase  shift  encountered  by  r.f. 

signal  In  propagating  to  target  and  back 
at  all  8  frequency  steps  and  various 
target  ranges  for  an  8-frequency  step 
frequency  radar. 


correspond  to  A,  Aexp(jTT/4),  Aexp(jTT/2),  ...  , 
Aexp(j7TT/4)  (See  Fig.  4.9).  The  real  components  of 
this  series  are  the  indicated  samples  of  the 
cosinusoidal  wave  shown  in  Fig.  4.10  at  R=Ra/8,  and  the 
imaginary  components  are  the  indicated  samples  of  the 
corresponding  sinusoidal  waveform  in  Fig.  4.11. 

Samples  of  waveforms  for  the  other  7  ranges  are  shown 
as  well  in  Figs.  4.10  and  4.11. 

Notice  that  the  frequency  of  the  sampled  waveform 
increases  as  the  range  increases  to  a  range  of  Ru/2. 

At  ranges  greater  than  Ru/2  the  frequency  decreases  and 
the  sign  of  the  sinusoidal  waveform  (i.e.,  the 
imaginary  component)  becomes  negative.  The  sign  change 
in  the  imaginary  component  makes  it  possible  to 
distinguish  between  ranges  of  Ru/2  <  z  <  Ru.  and 
0  <  z  <  Ru/2.  We  cannot  distinguish  ranges  greater 
than  Ru  from  ranges  less  than  Ru  due  to  the  sampling 
process . 

Assuming  that  the  return  amplitude  is  constant  at 
each  step,  the  return  for  a  point  target  at  a  range  of 
z  *  is 

hv=exp(-j0t)exp(  ( j4TTfc/c)z*  -hj{Znb/c)zs) 

N- 1 

x  E  Aexp(  (  j4TTnAf/c)z*  )  exp  (  -  j  2TT(  kn/N  )  ) 
n  =  0 

=  exp(-j0i)exp(  (  j4TTfc/c)Z*  +  jfeirfe. /c) Zs) 
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X 


Aexp(-j{  (N-U/2)  (2TT(2A£z*  /c-k/N)  ) 
x  sin(N/2(2TT(2Afz.*  /c-k/N)  )  )/sin(  1/2  (  2TT(  2££  /c-k/N)  )  ) 

=  Aexp(  20cZc+esZs 

x  sin(N/2(2TT(2Afzt/c-k/N)  )  ) /s  in  ( 1/2  (  2TI(  2^£zt /c-k/N )  )  ) 

(4.15) 

Notice  the  similarity  between  this  expression  and 
the  one  shown  for  PAIR;  only  the  2TTAf TZ  term  is  absent 
here.  As  was  the  case  with  PAIR,  amplitude  data  are 
easily  obtained.  Copolar  phase  is  determined  in  the 
same  manner  as  for  PAIR. 

The  sense  of  the  phasor  rotation  is  somewhat 
arbitrary  and  depends  on  the  phase  reference.  In  our 
radar  the  phasor  rotation  opposes  that  of  Iizuka  such 
chat 
N-l 

Hr,  =  Z  (skE0/z*  )exp(-j(  (4TI<fc+nAf  )/c)  (z0+fcAz)  ) 
k  =  0 

I f  we  set 

H'n=NHnexp(-j2(2n(fc  +n£>f)/c)  (z0  )  )  and 
h*  =  (SfcE0/ZK  )exp(  -  j  (  4TT/c  )  (  kfc  £z  )  ) ,  then 
N-l 

HV  =1  hkexp(  -  j  (  2TT(  kn/N  )  )  ) 
k=0 

N-l 

and  hk=(l/N)I  H  'n  exp  (  j  (  2TT(  kn/N )  )  )  .  (4.16) 

n  =  0 

That  is,  we  use  the  inverse  FFT  to  discriminate  target 
ranges  given  the  sampled  vector  measurements. 


4.3  Calibration  and  performance 

4.3.1  SOURCESCAT 

We  used  a  12"  Luneberg  lens  to  calibrate  this 
radar  during  the  1987  experiment  and  found  only  about 
±1  dB  variation  in  the  return  over  a  temperature  range 

Q  Q 

of  -25  C  to  5  C.  Delay  line  readings  were  taken  to 
track  the  variation. 

4.3.2  PAIR 

We  used  an  8"  Luneberg  lens  and  an  active  radar 
calibrator  (ARC)  to  obtain  calibration  readings  for 
this  radar  during  the  1988  experiment.  At  the  time  ve 
were  unaware  that  a  Luneberg  lens  does  not  work  well 
for  systems  with  small  bistatic  angles.  Thus  we  had  to 
rely  on  a  comparison  of  returns  from  the  lens  and  the 
not-vell-character ized  ARC  to  obtain  a  calibration 
value.  Variations  in  the  delay  line  readings  were 
observed  to  be  about  ±1  dB. 

During  the  1989  experiment  we  again  used  the  8" 
Luneberg  lens,  because  of  our  single  antenna 
configuration.  We  found  a  variation  of  about  ±2  dB  in 
the  return  over  a  temperature  range  of  -15°C  to  5°C. 
Delay  line  readings  were  taken  to  track  the  variation. 

The  phase  performance  for  this  radar  has  not  been 
well  characterized.  We  observed,  using  an  HP-3585 
spectrum  analyzer,  that  ve  were  able  to  adjust  the  PLL 
to  obtain  a  nearly  ideal  sine  function  return  for  a 


Phase  Difference  Between  HH  &  VV  Pol  for  Kulti-feor  ice 


Figure  4.12:  Uncalibrated  copolar  phase  measured  by 

PAIR  for  several  spots  on  a  smooth  saline 
Ice  surface. 


point  target.  We  also  found  that  the  radar  has 
reasonable  phase  stability;  figure  4.12  indicates  the 
copolar  phase  for  several  spots  on  a  smooth  ice  surface 
at  20°  incidence,  taken  over  a  one  hour  time  interval. 
Thus  ve  have  some  confidence  in  this  radar. 

The  performance  for  weaker  targets  and  long-term 
stability  have  not  been  determined. 

4.3.3  STEFRAD 

We  used  an  8"  Luneberg  lens  and  a  diplane 
consisting  of  2-10"xl0"  metal  plates  to  calibrate  this 
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radar  during  the  1989  experiment  and  found  about  ±1  dB 
variation  in  the  return  over  a  temperature  range  of 
-15°C  to  5°C.  Delay  line  readings  were  taken  to  track 
the  variation. 

We  have  been  able  to  study  the  phase  performance 
in  some  depth.  Figure  4.13  illustrates  measured 
backscatter  from  an  8"  sphere  at  C  band,  before  and 
after  subtraction  of  the  phase  response.  Figure  4.14 
shows  the  same  measurement  at  Ku  band.  Note  that  the 
system  response  is  reduced  by  nearly  40  dB  for  C  band. 
In  these  cases  the  sky  (system)  and  target  responses 
were  taken  within  a  minute  of  each  other,  however. 

The  C-band  radar  is  considerably  less  stable  than 
the  Ku-  band  radar  over  time.  With  the  Ku-band  radar 
we  were  able  to  remove  >30  dB  of  the  system  response  30 
minutes  to  1  hour  after  a  sky  reading  was  taken.  A  sky 
reading  for  C  band  was  good  for  less  than  5  minutes. 

A  simple  experiment  was  done  to  study  the  copolar 
phase  response:  VV  and  HH  polarization  responses  were 
determined  for  an  8"  metal  sphere.  A  known  phase  shift 
was  then  inserted  in  the  VV  channel  and  the 
measurements  were  repeated.  After  removing  the  system 
response  of  each  sphere  measurement  we  determined 

0VV~^|,K  —  ^vvmieri  —  0HH  ir>ier4 

^  ^Wnomtert  &  HH  *->«i«ser+  ^  (4.17) 

for  each  of  three  inserts.  Figure  4.15  shows  the 


Figure  4.14:  Return  from  an  8M  sphere  as  output  from 
Ku-band  STEFRAD  (a)  before  and  (b)  after 
removal  of  the  system  response. 
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Figure  4.15:  Results  of  C-band  copolar  phase 
•  experiment  using  inserts  of  known  phase 

shift. 
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Figure  4.16:  Results  of  Ku-band  copolar  phase 

experiment  using  inserts  of  known  phase 
shift. 
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results  for  c  band  and  Figure  4.16  the  results  for  Ku 

band  -  notice  that  we  plotted  the  phase  as  it  was 
determined  by  the  network  analyzer;  in  the  presence  of 
multiple  targets  it  is  necessary  to  determine  copolar 
phase  at  the  appropriate  peaks  of  the  IFFT. 

The  measurements  for  C  band  were  done  outdoors 
using  a  styrofoam  mount  for  the  sphere.  The  measured 
phase  was  generally  within  ±30°  of  the  predicted  phase. 
Errors  in  the  unwrapping  routine  are  responsible  for 
the  large  spikes;  notice  that  these  are  of  magnitude 
2TT.  Motion  of  the  mount,  due  to  a  slight  breeze, 
caused  most  of  the  remaining  variation. 

The  measurements  at  Ku  band  were  conducted  inside 
with  antennas  pointing  directly  into  a  large  metal  door 
which  was  partly  covered  with  microwave  absorbers.  The 
measured  phase  was  generally  within  ±20°  of  the 
predicted  phase  except  for  insert  1  (the  offset  for 
insert  1  can  be  attributed  to  a  faulty  connector). 
Errors  in  the  unwrapping  routine  are  responsible  for 
the  large  spikes,  the  same  as  for  the  C  band  data. 
Vibration  of  the  styrofoam  stand  by  the  room 
ventilation  and  interference  from  the  door  (the  door 
return  could  not  be  completely  removed  due  to  blockage 
from  the  sphere)  caused  most  of  the  remaining 


5.0  Analysis  of  ice  sheet  physical  parameters 

In  this  chapter  we  determine  parameters  to  be  used 
as  inputs  to  the  theoretical  models  discussed  in 
Chapter  2,  and  thus  prepare  the  groundwork  for  the 
comparisons  of  measurements  and  theory  presented  in 
Chapter  6.  The  ice  sheet  parameters  we  determined  are 
listed  in  Tables  3.2  and  3.3. 

5.1  Surface  roughness 

During  1987  we  estimated  rms  height  by  inspection 
and  used  contextual  information,  previously  reported 
surface  measurements,  and  fits  to  theory  to  estimate 
correlation  length.  These  techniques  were 
unsatisfactory  because  of  the  large  uncertainty. 

During  1988,  for  each  ice  sheet,  a  sample  was  cut 
from  the  ice  and  photographed.  The  resulting  surface 
profiles  for  smooth  and  rough  saline  ice  are  shown  in 
Figs.  5.1.  Although  good  profiles  were  obtained,  these 
results  were  still  lacking  in  two  respects:  (1)  only 
one  surface  profile  was  done  for  each  ice  sheet;  and 
(2)  the  block  was  removed  so  as  not  to  interfere  with 
subsequent  measurements;  i.e.,  a  block  was  cut  outside 
of  the  area  being  studied. 

During  1989,  a  wire  gauge  was  used  to  monitor 
roughness.  With  this  device  it  was  possible  to  obtain 
many  profiles  within  the  area  being  studied,  but, 
unfortunately,  with  reduced  accuracy.  Typical  surface 
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roughness  profiles  measured  with  this  device  are  shown 
in  Fig.  5.2a,  for  smooth  saline  ice,  and  in  Fig.  5.2b 
for  desalinated  ice. 

We  determined  rms  height  and  autocorrelations,  an 
the  corresponding  uncertainties,  for.  the  surfaces 
observed  during  1988  and  1989  .-  The  rms  height  we 
determined  as 

r  N  _  a-)®*5' 

<r*(l/N)E.  Izi  -z  ]  (  (5.1) 

^  i  =  1  J 

where  z^  are  surface  height  measurements  at  spatial 
intervals  of  about  3  mm  and 
N 

z=(l/N)I  Zj 
i  =  1 

The  autocorrelations  we  determined  as 
AM 

Rxx$  =  ( 1/N )  It  (  z^ -z  )  ( z.^  -z)] 
i~  1 

and  normalized  autocorrelations  as 
rxxi=Rxxi/Rxx0 

where  i  =  0,1, 2, . . . ,N-1  (5.2) 

As  a  first  approximation  we  fit  each 
autocorrelation  to  an  exponential  function  of  the  form 
exp(-x/l),  where  1  is  the  correlation  length  and 
1=0. 003k,  where  k  is  the  value  of  i  corresponding  to 
rxxi=l/e.  The  autocorrelations  and  corresponding  fits 
for  the  smooth  saline,  rough  saline,  and  desalinated 
ice  surfaces  of  1988  are  shown  in  Fig.  5.3  and  those 


Figure  5.2:  Surface  profiles  for  (a)  smooth  saline  ice 
sheet  and  (b)  desalinated  ice  blocks 
studied  during  the  CRREL'89  experiment. 
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Figure  5. 


3:  Computed  autocorrelations  and  miniraum- 
raean-3quare  exponential  fits  for  surfaces 
of  (a)  smooth  saline,  (b)  rough  saline, 
and  (c)  desalinated  ice  sheets  studied 
during  the  CRREL'88  experiment. 
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Figure  5.4:  Computed  autocorrelations  and  rainimum- 

raean-square-exponent lal  fits  for  surfaces 
of  (a)  smooth  saline  ice  sheet,  and  (b) 
desalinated  ice  blocks  studied  during  the 
CRREL'89  experiment. 


1 


for  the  smooth  saline  and  desalinated  ice  surfaces  of 
1989  are  shown  in  Fig.  5.4.  The  correlation  lengths 
for  each  surface  are  listed  in  Tables  3.2  and  3.3. 

We  then  sought  to  determine  standard  deviations 
for  the  rms  heights  and  correlation  lengths,  which  are 
nonzero  because  of  finite  sample  sizes,  so  that  bounds 
for  theoretical  predictions,  accounting  for  uncertainty 
(i.e,  rms  height  ±1  standard  deviation  and  correlation 
±1  standard  deviation)  in  our  parameter  estimates, 
could  be  obtained.  Since  the  surface  profiles  of  1988 
consisted  of  1  profile  for  each  surface  it  was  not 
possible  to  obtain  those  statistics  from  the  data. 

Hence  we  performed  Monte  Carlo  simulations  to  obtain 
more  profiles,  making  the  following  assumptions  about 
the  statistical  distributions  of  the  measured  surface 
heights.  We  assumed  that  the  surface  height  data  were 
wide-sense  stationary,  gaussian  distributed,  and 
exponentially  correlated,  and  plotted  histograms  of  ti.e 
measured  data  to  verify  that  these  assumptions  were 
reasonable  . 

To  generate  samples  with  these  statistical 
properties  one  can  pass  samples  from  a  white  gaussian 
distribution  (i.e.,  w(x))  through  a  linear  first-order 
system  of  the  form 

dz(x)/dx  +  Az  ( x )  =  Bw( x )  (5.3) 

where  A  =  inverse  of  the  correlation  length 


B  is  adjusted  so  that  B*  =  rms  surface  height 

**  =  E{ w2 ( x ) } 

z(x)  is  the  output. 

The  statistical  parameters  of  z(x)  are  related  to  the 
system  parameters.  A,  B  and  v(x)  by 
E{z2(x) } =B2  a2 
and 

E{z(x)z(x+x' ) } =  B2 a2 exp( -A  lx  1  I ) . 

For  convenience  ve  implemented  a  discrete  version  of 
this  system,  i.e., 

zk+1  =Fzk+Gvk  (5.4) 

on  a  personal  computer.  The  continuous  and  discrete 
systems  give  approximately  the  same  results  for  small 
spatial  sampling  interval,  Ax,  if 
F=exp ( - AAx ) , 

G= ( 1/A )( 1-exp ( -AAx )) B,  and 

vk=v(x)/yST  (5.5) 

[ Brown,  198  3  ]  . 

For  each  surface  of  1988  we  simulated  20  profiles 
the  length  of  the  measured  profile  to  obtain  estimates 
of  the  standard  deviations  of  the  rms  height  and 
autocorrelation.  Because  the  autocorrelations  of  10 
measured  profiles  were  averaged  to  obtain  an 
autocorre lat ion  for  each  surface  of  1989,  we  simulated 
200  profiles  the  length  of  the  measured  profiles, 
computed  the  autocorrelations  for  each  profile. 
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averaged  them  by  10's,  and  used  the  20  averages  to 

obtain  the  standard  deviation  limits  of  the 
autocorrelation  computed  from  measurements.  Standard 
deviation  of  the  rms  height  was  similarly  computed  from 
20  averages  of  10  rms  heights. 

The  autocorrelations  and  estimated  ±1  s.d.  limits 
for  each  surface  are  shown  in  Fig.  5.5  for  the  1988 
surfaces  and  in  Fig.  5.6  for  the  1989  surfaces.  Only 
the  rough  ice  surface  autocorrelation  (Fig.  5.5b) 
deviates  greatly  from  its  limits,  implying  that  it  may 
not  be-  exponential.  As  a  first  approximation  the 
exponential  form  is  considered  because  there  are  not 
enough  data  points  to  determine  if  the  secondary  peaks 
of  Fig.  5.5b  are  statistically  significant.  Notice 
that  the  uncertainties  for  1989  tended  to  be  smaller 
than  those  for  1988  because  more  profiles  were  obtained 
for  each  surface  in  1989. 

5..  2  Internal  structure 

Thin  section  photographs  of  saline  ice  have  been 
supplied  by  CRREL .  In  addition,  during  1989,  we 
photographed  small  samples  and  horizontal  and  vertical 
thick  sections  from  near  the  top,  middle  and  bottom  of 
the  desalinated  ice  blocks  (See  Chapter  3). 

From  the  thick  section  photographs  ve  located  as 
many  of  the  air  bubbles  as  we  could  identify  and 
highlighted  them  with  a  marking  pen.  We  then  made 
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Figure  5.5:  Computed  autocorrelations  and  ±  1  standard 
deviation  limits  obtained  from  Monte  Carlo 
simulations  for  surfaces  of  (a)  smooth 
saline,  (b)  rough  saline,  and  (c) 
desalinated  ice  sheets  studied  during  the 
CRREL’88  experiment. 
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Figure  5.6:  Computed  autocorrelations  and  £  1  standard 
deviation  limits  obtained  from  Monte  Carlo 
simulations  for  surfaces  of  (a)  smooth 
saline  ice  sheet,  and  (b)  desalinated  ice 
blocks  studied  during  the  CRREL'89 
experiment. 


Figure  5.7 j  Bubbles  In  a  3— mm  thick  section  taken  from 
the  top  of  a  sample  of  a  desalinated  Ice 
block  studied  during  the  CRREL’89 
experiment,  as  visually  Identified  In  a 
thick-section  photograph. 


photocopies  of  the  photographs,  used  vhite-out  where 
necessary  to  eliminate  the  background  and  made  more 
photocopies  (See  Fig.  5.7).  We  then  sent  these  to 
CRREL  where  an  image  processor  was  used  to  determine 
correlation  lengths  of  the  spacings  between  air 
bubbles.  The  analysis  of  each  processed  thick  section 
was  performed  assuming  air  bubbles  were  of  value  1 
(i.e.,  the  relative  dielectric  constant  of  air)  and 
the  background  was  of  value  3  (i.e.,  the  relative 


dielectric  constant  of  ice).  Mean  and  standard 


Upper  3 

mm 

horiz  sect 

2.00 

1.77 

Middle  3 

mm 

horiz  sect 

2.05 

2.65 

bottom  3 

mm 

horiz  sect 

2.17 

2.17 

Upper  3 

cm 

vertical 

2.00 

1.37 

Middle  2 

cm 

vertical 

2.79 

1.48 

Bottom  2 

cm 

vertical 

4.28 

2.11 

LI  and  L2  are  correlation  lengths  in  mm  on  major  and 
minor  axis  of  correlation  ellipse. 


Table  5.1:  Correlation  lengths  of  air  bubbles  in  a 
sample  taken  from  a  desalinate  ice  block 
studied  during  the  CRREL'89  experiment, 
computed  from  thick-section  photographs  that 
were  processed  as  shown  in  Fig.  5.7. 


deviation  of  the  image  and  a  corresponding  subimage, 
near  the  center  of  the  image,  were  compared  -to  be  sure 
that  the  statistics  were  the  same.  Then  the  subimage 
was  moved  across  the  image  to  obtain  a  two-dimensional 
autocorrelation  such  as  that  shown  in  Fig.  5.8a  for 
the  upper  3  mm  horizontal  section  and  in  Fig.  5.8b  for 
the  middle  2  cm  vertical  section  (Autocorrelations  of 
the  other  sections  are  similar  to  one  of  these  two). 

The  correlation  lengths  (i.e.,  spacing  between 
points  along  the  plane  of  the  thick  section  which 
results  in  a  normalized  autocorrelation  of  1/e)  along 


Figure  5.8:  Autocorrelation  of  air  bubbles  in  (a)  a  3- 
mm  thick  horizontal  section  of  the  top  and 
(b)  3-mm  thick  vertical  section  near  the 
middle  of  a  sample  taken  from  the 
desalinated  ice  blocks  studied  during  the 
CRREL'89  experiment,  computed  from  thick- 
section  photographs  that  were  processed  as 
Shown  in  Fig.  5 . 7  .  (Courtesy  of  CRJtEL) 


the  major  and  minor  axes  of  autocorrelation  plot  were 
then  determined  and  are  listed  in  Table  5.1.  Notice 
that  the  correlation  lengths  of  the  horizontal  sections 
are  nearly  symmetric  whereas  those  of  the  vertical 
sections  have  correlation  length  ratios  (major 
axis/minor  axis)  of  1.5  -  2,  increasing  with  depth  -  a 
number  of  vertically  oriented  ellipsoidal  bubbles  in 
the  lower  half  of  the  ice  were  responsible  for  this. 

From  the  plots  in  Fig.  5.8  it  seems  reasonable  to 
estimate  the  autocorrelation  function  as  exponential. 
Simulations  similar  to  those  done  for  the  surfaces 
would  be  informative  but  were  not  done. 

A  volume  fraction  profile  of  the  air  in  the 
desalinated  ice  was  estimated  using  measured  bulk 
density/  air  bubble  size  distributions,  and  approximate 
bubble  densities  (number  of  bubbles/cm1).  We  used  the 
air  bubble  distributions  and  bubble  densities,  computed 
from  processed  horizontal  thick  sections  taken  from 
near  the  top,  middle  and  bottom  of  the  desalinated  ice, 
to  determine  a  relative  volume  fraction  profile.  We 
then  used  these  values  and  the  bulk  density  measurement 
of  0.844  g/cm3,  assuming  the  ice  consisted  only  of  pure 
ice  and  air,  to  determine  absolute  values  for  the 
volume  fraction  profile.  The  profile  listed  in  Table 
5.2  was  obtained  assuming  that  the  ice  consisted  of 
three  equal-size  layers,  with  bubble  distribution  and 
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Case 


Volume  fraction  of  air 


Top  4 . 5  cm  20.3% 

Middle  4.5  cm  2.2% 

Bottom  4.5  cm  0.12% 


Table  5.2:  Volume  fraction  profile  of  air  in  the 

desalinated  ice  blocks  observed  during  the 
CRREL'89  experiment,  computed  from  estimated 
air  bubble  density  and  size  distributions, 
and  a  bulk  density  measurement. 


density  in  each  layer  being  the  same  as  that  measured 
for  the  corresponding  horizontal  section. 

5.3  Snow  cover 

One  useable  snowfall  occurred  during  the  1987 
experiment.  Unfortunately  the  snow  parameters  were  not 
measured.  One  snowfall  also  occurred  during  1989.  The 
parameters  of  this  snow  are  listed  in  Table  3.2. 

The  surface  profile  of  the  1989  snow  was 
estimated  from  a  photograph.  A  plot  of  roughness  and 
correlation  function  for  the  snow  cover  of  1989  are 


shown  in  Figs.  5.9. 


6.0  Data  Analysis  and  comparison  with  theory 

We  used  an  analog  spectrum  analyzer  to  convert  the 
SOURCESCAT  output  signal  into  a  more  useable  range- 
discriminated  form  (frequency  domain).  For  PAIR  and 
STEFRAD  a  Discrete  Fourier  Transform  (DFT)  program 
performed  the  same  function.  Note  that  the  output  of 
STEFRAD  is  equivalent  to  the  time  domain  outputs  of 
SOURCESCAT  and  PAIR  -  see  Chapter  4  for  a  discussion  of 
point  target  returns  for  each  of  the  radars. 

The  BASIC  routines  ve  used  to  study  and  process 
the  SOURCESCAT  data  were  slow  and  relatively 
inflexible,  so  ve  wrote  a  comprehensive  data  analysis 
program  in  TURBO  Pascal  to  study  and  process  PAIR  and 
STEFRAD  data.  The  most  important  components  of  our 
data  analysis  software  are:  (1)  a  flexible  CRT  graphics 
display  routine;  (2)  an  efficient  FFT  algorithm;  and 
(3)  an  accurate  integration  routine.  Other  important 
components  of  our  software  include  routines  to  subtract 
system  response,  unwrap  phase,  and  determine  phase 
differences.  This  package  enabled  us  to  study 
amplitude  and  phase  of  the  return  signal,  in  both  the 
time  and  frequency  domains,  and  to  integrate  over 
selected  intervals. 

6.1  Computation  of  <r° 

In  determining  <r° ,  it  is  necessary  to  compute  the 
total  power  returned  from  the  target.  Total  target 


power  Is  obtained  by  integrating  over  the  appropriate 
segment  of  the  range-discriminated  signal  (i.e., 
frequency  domain).  (Note  that  range  discrimination  is 
necessary  even  in  the  apparent  absence  of  other  targets 
due  to  radar  system  reflections  and  feedthroughs  as 
shown  in  Figs.  6.1).  The  integration  for  SOURCESCAT 
was  performed  as 
n 

P  =  I  P ( «| )  (6.1) 

i=m 

where  P(«i)  =  power  density  at  frequency  “i  as  obtained 

from  the  analog  spectrum  analyzer 
=  low  frequency  limit  of  target  return 
=  high  frequency  limit  of  target  return 
The  integration  for  PAIR  and  STEFRAD  was  performed  as 
n  n 

P  =  I  X(«x)X*(« i)  =  I  I X ( ) 12  (6.2) 

i=m  l=m 

where  X(«*»)  is  the  DFT  of  the  radar  output  and  and 
are  defined  above. 

Once  total  power  values  were  obtained,  standard 
narrow-beam  techniques  were  used  to  obtain  a0;  i.e., 

Prt  »CPtGt  Aecr°A/(  4TT)2  R*  (6.3) 

where  P,.t  =  total  return  power  from  target 
P*  =  transmitted  power 
Gt  =  antenna  gain 
Ae  =  antenna  effective  aperture 

=  characteristic  cross-section  of  target 
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Figure  6.1:  Returns  from  target  and  system  as  seen  in 
range-discriminated  outputs  of  (a)  X-band 
SOURCESCAT  (b)  Ku-band  STEFRAD. 
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A  =  Illuminated  area  of  target 
R  =  slant  range  to  the  target 

C  =  constant  to  account  for  system  parameters 

not  in  radar  equation;  e.g.,  system 
transfer  functions  and  loss 
Since  many  of  these  parameters  were  not  accurately 
known,  measurements  of  an  internal  reference  (delay 
line)  were  used  to  monitor  system  changes  to  allow 
comparing  target  returns  with  return  from  a  target  of 
known  cross-section,  a.  Let 
P0Lt  =K2K^L  and 

Prt  =  Kt  Kfc  a°  A/R*  ,  (6.4) 

assuming  does  not  change  from  target  measurement  to 
delay  line  reading,  and 
where  Prt  is  defined  above  and 

P0L*  =  delay  line  return  power  at  the  time  of 
target  measurement 

Kx,Ka  =  transfer  functions  of  stable  system  paths 
of  radar  in  radar  mode  and  delay-line 
mode  respectively 

K-t  =  transfer  function  of  time-  and  temperature- 
varying  system  path  at  time  of  target 
measurement 
L  =  delay  line  loss 
Rt  =  slant  range  to  the  target 


PDL.  =K2kcL  and 

P~  K  c  a/R*  ,  (6.5) 

assuming  Kc  does  not  change  from  calibration 
measurement  to  delay  line  reading.  Ka,K2,  and  L  are 
defined  above  and 

P0Lc  =  delay  line  power  at  the  time  of 
calibration  measurement 
Pj-C  =  power  return  from  known  target 
Kc  =  transfer  function  of  time-  and  temperature- 
varying  system  path  at  time  of  calibration 
measurement 

R<.  =  slant  range  to  calibration  target 
o  =  cross-section  of  known  target, 
then  a0  can  be  related  to  a  via  a  series  of  ratios. 
prt  /Pre  =<K*/Kc  )  (o°A/o)  (R*  /R*  ) 

=  (PDLt  /PBtc  )(*°A/a)(Rj  /R*  ) 


and 


•  °G  =  (Prt/PrcHPDU  /PDLl  )(a/A)(R*t  /R’t  )  (6.6) 

Once  narrow-beam  values  were  obtained,  a  routine 
by  Kim  (1982]  was  used  to  determine  if  a°  variation 
within  the  beam  was  sufficiently  small  so  that  the 
narrow-beam  approximation  was  valid.  In  addition,  a 
routine  by  Wang  (1989]  was  used  to  make  necessary 
corrections  for  antenna  separation. 

£  6.1.1  Consideration  of  errors 

The  calibration  targets  used  were  8"  and  12” 
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Luneberg  lenses,  a  diplane  constructed  of  2  10"xl0" 
metal  plates,  and  an  active  radar  calibrator.  Returns 
from  the  diplane  (a=4TTA*ff  /x2=  16.2  dB  at  C  band  and 
23.3  dB  at  Ku  band)  were  used  to  verify  the  predicted 
cross-section  of  the  8"  Luneberg  lens  (1  to  2  dB  less 
than  4TI*r  / \  ,  or  6  dB  at  C-band  and  12  dB  at  Ku  band). 
The  predicted  cross-section  of  the  12"  Luneberg  lens  (1 
or  2  dB  less  than  4TT3r*/x2',  or  17  dB)  was  not  properly 
verified  -  this  i3  not  crucial  to  the  discussion  to 
follow  in  that  only  a  qualitative  description  of  the 
data  depending  on  calibration  with  this  lens  is  given 
(i.e.,  the  CRREL'87  data). 

The  calibration  measurements  of  the  1988-1989 


experiment  were  done  with  the  8"  Luneberg  lens  and  are 
accurate  to  within  about  ±1  dB,  assuming  a  calibration 
range  uncertainty  of  ±10  cm  and  allowing  for  a  few 
tenths  of  a  dB  uncertainty  to  account  for  pointing 
error  between  the  antenna  and  lens.  The  calibration 


measurements  of  the  1986-1987  experiment  were  done  with 
the  12"  lens  and  are,  somewhat  arbitrarily,  assigned  an 
uncertainty  of  ±2  dB  to  allow  for  uncertainty  in  the 
12"  lens  cross-section.  The  C-band  calibrations  of  the 
1987-1988  experiment  were  done  with  the  8"  Luneberg 
lens  and  have  a  somewhat  wider  uncertainty  than  those 
of  the  other  2  experiments  (i.e.,  ±3  dB)  because  ve 
later  discovered  that  a  Luneberg  lens  is  not  a  suitable 
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calibration  target  for  a  system  with  a  bistatic  angle 

between  transmit  and  receive  antennas  (See  Pig.  6.2). 
The  lens  beamvidth  at  C  band  is  about  twice  that  shown 
in  Fig.  6.2  since  the  measurements  shown  in  the  curve 
were  done  at  X  band;  thus  our  calibration  readings 
would  have  been  about  5  dB  below  the  peak  return  of  the 
lens.  The  uncertainty  was  estimated  to  be  £3  dB  to 
allow  for  uncertainty  in  the  bistatic  angle  -  notice 
the  steepness  of  the  curve  in  the  vicinity  of  -5  dB. 

Other  uncertainties,  including  target  range  and 
error  due  to  use  of  the  narrow-beam  approximation, 
account  for  a  maximum  error  of  ±1  dB.  These 
uncertainties  are  referred  to  as  target  measurement 
uncertainty  in  the  following  discussion. 

Another  consideration  in  estimating  o°  is 
uncertainty  due  to  fading  -  because  the  number  of  data 
points  obtained  in  each  case  is  too  small  to  be  used  to 
estimate  the  actual  distribution,  a  typical  exponential 
distribution  is  assumed.  During  the  1987  campaign  we 
obtained  a  minimum  of  10  uncorrelated  samples  for  each 
measurement.  Assuming  Rayleigh  fading  l±i*/a  =  1.19, 
-1.65  dB.  During  the  1988  and  1989  campaigns  we  were 
able  to  obtain  only  5  uncorrelated  samples,  in  most 
instances,  and  3  in  some  cases.  Assuming  an  exponential 
distribution,  l±i*/o  =  1.6,  -2.7  dB  for  5  uncorrelated 
samples  and  =  2.0,  -3.74  dB  for  3  uncorrelated  samples. 
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To  consider  the  combined  effect  of  calibration  and 
target  measurement  uncertainties  and  fading  on  ®°,  we 
express  equation  (6.6)  In  the  form 

=  kFPrt  ®  (6.7) 

where  k  =  constant 

F  =  factor  to  account  for  target  measurement 
uncertainty 

Pr*  =  return  power  obtained  from  a  average  of  a 
designated  number  of  independent  samples 
a  =  lens  cross  section. 

We  then  make  the  assumption  that  F,  Prt  and  a  are 
statistical  parameters  with  F  accounting  for  target 
measurement  uncertainty,  Pc*  accounting  for  fading,  and 
a  accounting  for  calibration  uncertainty.  The  combined 
uncertainty  can  then  be  estimated  by  expanding  101og(®°) 
in  a  Taylor  series  around  mf,  MPrt  and  m„.  and  taking 
expected  values.  From  Breipohl  (19701: 

If  Y  =  g(X,,  X2,  ...,  X„)  and  Mi  =  EIX*]  then 


Y  ~  g(Mj  ,m2. 


and 


n 


.,M«)  +  I  t£2.(Mi,Mi, 

i=l  tfX i 


. , M*  )  1  ( X j  ~Mi ) 


n 


o*  ~  1  [  £ol(m,  ,m2  , .  . .  ,Mn)  ]  <Ti  (F,Prt.,®  independent) 

1=1  6X i 


(6.8) 


So  i  f 
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* 


101og(o°)  =  101og(k)+101og(F)+101og(Prt )+101og(*) 
then 

a*  =  (101og(e)  )*  I  +  (o^  /it^  >*) 

(in  dB) 

^  The  estimated  total  uncertainty  values  were 

obtained  as  ±o  and  are  listed  in  Table  6.1.  They  are 
also  included  on  plots  containing  measured  data  points 
•  as  vertical  bars. 

uncertainty 

target  Fading  Total 


cal . 

1987 

• 

all  surfaces 

±2 

±1 

-1.7, +1.2 

±3 

1988 
smooth 
saline  ice 

±3 

±1 

-2. 7, 1.6 

±4.9 

• 

rough  saline 

±3 

±1 

-1.7,1. 2 

±4.7 

ice 

desalinated 

ice 

±3 

±1 

-2. 7, 1.6 

±4.9 

• 

1989 

smooth  and 
3nov  covered 
smooth  saline 

±1 

±1 

-2. 7, 1.6 

±2.5 

ice 

desalinated 

±1 

±1 

-3.7,2 

±3 

• 

ice  blocks 

Table  6.1:  Estimated  measurement  and  fading  uncertainty 
(in  dB)  for  the  three  years  of  experiments. 

6.2  Results 

In  this  section,  the  processed  data  are  presented. 
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In  particular  we  treat  the  effects  of  surface 
roughness,  snow  '•over,  and  volume  inhomogeneities.  The 
Kirchhoff  scalar  approximation  and  Small  Perturbation 
Method  surface  scatter  models  and  the  Effective  Medium 
volume  scatter  theory  (see  Chapter  2)  are  used  in  the 
comparisons  between  results  and  theoretical 
predictions . 

6.2.1  Smooth  saline  ice 

Backscatter  from  smooth  saline  ice  at  X  band  (on 
Jan.  17,  1987)  is  illustrated  in  Fig.  6.3.  Notice  that 

and  a°H  are  not  statistically  different.  We 
determined  the  penetration  depth  of  the  ice  to 
determine  if  measurable  volume  backscatter  was  possible 
and  used  the  Kirchhoff  surface-scatter  model,  with 
estimated  surface  parameters  as  inputs,  to  determine  if 
surface  scatter  alone  could  account  for  our  results. 

We  determined  the  penetration  depth  from  estimates 
of  the  ice  dielectric  constant,  using  the  model  of 
Frankenstein  and  Garner  [19671  with  the  salinity 
profile  of  Fig.  3.5,  and  assumed  linear  temperature 
variation  as  a  function  of  ice  depth  (-15°C,  or  air 
temperature,  at  the  surface  to  near  0°C  at  the  water 
surface).  Returns  from  the  ice-water  interface  (at  6 
cm  depth)  were  negligible,  since  they  were  attenuated 
by  more  than  12  dB.  Other  contributions  from  within 
the  ice  were  not  ruled  out;  e.g.  return  from  the 
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•  Smooth  Ice  -  V  Pol  6 
x  Smooth  Ice  -  H  Pol  I 


0  10  20  30  40  50 

Angie  of  Incidence 

Figure  6.3:  Backscatter  cross-section  of  smooth  (l.e., 
ms  height  <0.1  cm)  saline  ice  sheet  at  X 
band,  as  determined  from  data  obtained 
during  the  CRREL'87  experiment. 


interface  between  the  more  granular  and  less  granular 
ice  at  a  depth  of  about  3  cm  (see  Fig.  3.4). 

We  used  the  estimated  rms  surface  height  and 
adjusted  values  of  the  correlation  length  in  the 
Klrchhoff  surface  scatter  model.  By  adjusting  the 
correlation  length  (to  a  reasonable  value)  we  were  able 
to  obtain  predicted  values  that  were  within  a  few  dB  of 
the  measured  data. 

The  range-discriminated  radar  output  data  from 
which  Fig.  6.3  is  derived  showed  that  the  returns  from 
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the  surface  were  broader  than  would  be  expected  for 
surface  scatter  only.  We  performed  a  simulation  to 
show  that  this  broadening  may  have  been  due  to  a  small 
volume  scatter  component.  The  point  target  return  is 
contrasted  to  the  ice  surface  return  in  Pig.  6.4a, 
where  the  footprint  on  the  ice  surface  was 
approximately  a  circle  of  12"  radius.  To  simulate  the 
broadening  effect  (See  Fig.  6.4b)  we  assumed  that  a°  at 
50°  Incidence  was  due  solely  to  volume  scatter,  that 
the  narrow-beam  approximation  was  valid,  and  that  the 
volume  scattering  coefficient,  av,  was  constant  such 
that  (from  Ulaby  et  al.,  (19821) 

m 

a°(50°)=a°  (50°)  =  (T4(50°  )  j*«vexp(  -2«<z  )dz  lcos0'  (6.9) 

0 

where  *  =  power  absorption  coefficient  of  the  ice 
calculated  using  the  estimate  dielectric 
constant 

z  =  distance  of  penetration  into  the  ice 

T  =  transmission  coefficient  at  air-ice  interface 
assuming  a  smooth  interface 

G'  =  angle  of  refraction  in  saline  ice 
=  volume  scatter  contribtion  to  ao 
We  then  divided  the  simulated  target  into  range  cells, 
determined  the  surface  and  volume  contributions  from 
each  range  cell  to  a0 (30°)  =  (30°)  +  (30°),  using 

®°( 30° )  obtained  from  measurements  and  equation  (6.9) 
to  determine  a  °  .  This  information,  the  known  slant 


Figure  6.4:  (a)  measured  backscatter  from  artificially 
grown  sea  ice  at  X  band  and  30  degrees 
Incidence,  contrasted  with  point  target 
return  (b)  simulation  of  relative 
contributions  of  surface  and  volume 
scatter  to  backscatter  from  artificially 
grown  sea  ice  at  X  band  and  30  degrees 
incidence. 


range.  Illuminated  area  and  antenna  beamvidth  were  used 

int  the  radar  equation  to  calculate  the  return  power 
as  a  function  of  range. 

We  concluded  that  the  return  from  this  ice  surface 
was  dominated  by  surface  scatter  but  that  a  small 
volume-scattering  component  was  probably  present. 

The  ice  structure  and  dielectric  profile  on  Jan. 
21,  1987  were  nearly  the  same  as  they  were  on  Jan.  17. 
The  major  differences  in  the  ice  sheet  were  a  6  cm 
increase  in  thickness  and  a  modified  roughness.  The 
comments  concerning  volume  scatter  made  for  the  3heet 
of  Jan.  17  are  applicable  here.  We  again  assumed  that 
surface  scatter  dominated  the  return  due  to  the 
relatively  high  absorption  of  the  ice.  We  used  the 
Small  Perturbation  Model  with  the  estimated  rms  height 
and  a  correlation  length  adjusted  to  achieve  a  best  fit 
with  the  data  (see  Fig.  6.5).  Notice  how  the  slope  of 
the  a° ( 0 )  curves,  especially  the  HH  curve,  differed 
from  those  of  Jan.  17  (Fig.  6.3). 

According  to  the  Small  Perturbation  Method,  the 
change  in  backscatter  as  a  function  of  incidence  angle 
that  occurred  from  Jan.  17  to  Jan.  21  can  be  explained 
by  a  change  in  correlation  length  (Figure  6.6a 
illustrates  the  effect)  and  a  slight  change  in  rms 
height  (See  Fig.  6.6b),  or  by  a  change  in  correlation 
function.  The  latter  possibility  cannot  be  ruled  out 
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Figure  6.5:  Backscatter  cross  -section  of  textured 

saline  ice  sheet  at  X  band/  as  determined 
from  data  obtained  during  the  CRREL'87 
experiment. 


since  we  textured  the  Ice.  The  observed  difference 
between  VV  and  HH  polarizations  in  Fig.  6.5  is 
predicted  by  the  Small  Perturbation  Method.  The  change 
in  this  difference  between  Jan.  17  and  Jan.  21  cannot 
be.  With  the  first-order  Small  Perturbation  Method, 
only  a  change  in  dielectric  constant  can  change  the 
difference  between  VV  and  HH,  and  the  effect  is  small 
(tenths  of  a  dB  as  €'  changes  from  2.5  to  4). 

Backscatter  from  a  smooth  saline  ice  sheet  at  C 
band  (Feb.  3,  1989)  is  shown  in  Fig.  6.7.  The  2 
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Figure  6.7:  Backscatter  cross-section  of  smooth  saline 
ice  sheet  at  C  band*  as  determined  from 
data  obtained  during  the  CRREL'89 
experiment,  compared  with  Small 
Perturbation  Method  predictions. 


dotted  curves  indicate  the  maximum  and  minimum  Small 
Perturbation  predictions  of  oo(0)  based  on  the  rms 
height  ±  1  standard  deviation  and  the  correlation 
length  ±  1  standard  deviation  of  the  surface.  The 
dashed  line  indicates  the  noise  level  of  the  STEFRAD 
radar.  Notice  that  this  radar  is  incapable  of 
measuring  the  predicted  return. 

The  discrepancy  between  measurement  and  theory  can 
be  attributed  to  volume  scattering,  radar  system 
effects,  and  uncertainty  in  surface  characterization. 
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To  illustrate  the  volume  scattering  effect. 


consider  the  volume  scattering  component  as  an 

equivalent  rough  surface  -  this  has  some  merit  in  light 

of  the  abrupt  change  in  dielectric  constant  at  depths 

of  2-4  cm  observed  for  Feb.  3,  1989  (see  Fig.  3.20a). 

Assuming  the  proposed  model  «r°  is  computed  as 

«e=Ta(«:q(iiv,ul.£ace  *exp(-2etdsece' )  )  +  «i\e  (6.10) 

where  «  =  power  absorption  coefficient  of  the  ice 

calculated  using  the  estimated  dielectric 
constant 

T  =  the  air-ice  transmission  coefficient 
assuming  a  smooth  interface 

a  =  the  power  absorption  coefficient  determined 
from  the  estimate  dielectric  constant 

d  =  the  depth  of  the  equivalent  surface  in  the 
ice 

S'  =  angle  of  refraction  in  saline  ice 

°ice  =  radar  cross-section  of  the  ice  surface 
determined  from  the  Small  Perturbation 
Method 


o“,ulv  surface  =  radar  cross-section  of  the 

equivalent  surface  determined  f 
the  Small  Perturbation  Method 

The  application  of  this  model  to  the  saline  ice 
sheet  of  Fig.  6.7  is  illustrated  in  Fig.  6.8,  with 
equivalent  rough  surface  parameters  of  z 

/zlrte*  i  =1.2,  ormJ  =0.5  cm,  and  1  =  1.8  cm.  Notice  that 
this  model  can  account  for  the  level  of  most  of  the 


rom 


data  points. 

Radar  system  effects  are  shown  in  Fig.  6.7  and 


Figure  6.8:  Backacatter  cross-section  of  smooth  saline 
ice  sheet  at  C  band,  as  determined  from 
data  obtained  during  the  CRREL'89 
experiment,  compared  with  predictions  of  a 
simple  layer  model. 


include  system  noise  and  coherent  contributions  due  to 
the  wide  beamvidth  of  the  antennas.  As  can  be 
observed,  these  effects  are  still  several  dB  below  the 
measured  values. 

Uncertainties  due  to  the  statistical  nature  of  the 
ice  surface  were  discussed  in  Chapter  5.  Uncertainty 
in  surface  roughness  measurements  contributes  a  less 
well-character izeable  uncertainty  that  can  be  large  for 
smooth  surfaces.  For  example,  from  the  photograph  we 
used  to  determine  roughness  for  one  of  the  1988  ice 
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sheets,  another  group  determined  a  »rm4  value  15% 
higher  than  our  armS  +  1  standard  deviation.  But  even 
their  estimates  give  theoretical  predictions  that  are 
about  6  dB  lover  than  vhat  ve  measured. 

We  conclude  that  volume  scatter  contributes 
measurably  to  backscatter  from  smooth  saline  ice  at  C 
band.  We  have  not  been  able  to  quantify  the  effect 
however . 

6.2.2  Snow  cover  effects 

We  noticed  that  a  deep  layer  of  snow  (15  cm) 
caused  considerable  volume  scatter  at  X  band;  an 
isolated  experiment  was  done  to  illustrate  this.  With 
no  snow,  a  dominant  peak  corresponding  to  the  ice 
surface  was  observed  (Fig.  6.9a).  After  15  cm  of  loose 
snow  was  piled  on  the  ice  surface  a  relatively  broad 
spectral  return  was  observed  (Fig.  6.9b).  This  was 
indicative  of  volume  scatter  from  the  snow  particles  in 
the  air-snov  mixture  (the  snow  was  dry  and  the 
temperature  was  about  -9°C).  Then  the  snow  was  packed 
to  a  depth  of  7.5  cm  and  more  measurements  were  done. 
The  packing  increased  the  homogeneity  of  the  layer,  as 
evidenced  by  the  decrease  in  spectral  broadening; 
because  of  the  decrease  in  volume  scatter  peaks 
corresponding  to  the  snow  and  ice  surfaces  became 
noticeable  (Fig.  6.9c). 

At  X  band  6.5  cm  of  dry  snow  had  little  effect  on 


(taga  (a) 


Figure  6.9:  Effect  of  snow  cover  as  seen  in  range- 

discriminated  output  of  X-band  SOURCESCAT 
during  the  CRREL'87  experiment.  (a) 
return  from  bare  smooth  saline  ice  surface 

-  peak  at  slant  range  of  about  3.6m 
corresponds  to  surface,  (b)  broadening 
effect  caused  by  volume  scatter  within  15- 
cm  snow  cover,  and  (c)  return  after  15  cm 
of  loose  snow  packed  to  thickness  of  7.5cm 

-  peak  at  3.4m  corresponds  to  snow  surface 
and  peak  at  3.8m  corresponds  to  ice 
surface . 
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as  was  determined  during  the  1987  campaign.  No 
significant  contribution  from  the  snow  cover  was 
evident  at  either  polarization  (See  Figs.  6.10).  The 
properties  of  the  snow  covering  the  surface  are  not 
known,  except  that  the  snow  was  dry.  Assuming  that  the 
snow  was  of  density  0.2  g/cu  cm  (i.e.,  /»s  =  0.2),  and 
using  the  empirically-determined  expressions  of  Ulaby 
et  al  [1982],  which  relate  extinction  coefficient  and 
volume  scattering  coefficient  to  snow  density  at  9  GHz 
(i.e.,  k«=0.77/»s  and  av=0.4/»s),  the  volume  scatter 
contribution  of  the  snow  is  estimated  from 

=(  ovcose/2k<?  )  ( 1-1/L*  (  0)  )  (6.11) 

where  ke  =  extinction  coefficient 

L(6)  =  loss  of  the  snow  layer  at  incidence  angle  0 
ke ,  <r„  are  defined  above 

to  be  approximately  -24  dB.  Since  this  is  6  to  8  dB 
lower  than  bare  ice  a(0),  it  is  not  surprising  that  no 
significant  contributi  -n  from  the  snow  was  evident. 

That  is,  the  ice  surface  was  sufficiently  rough  (see 
parameters  in  Fig.  6.3  )  that  surface  scatter  dominated 
over  volume  scatter  from  the  snow. 

During  1989  we  studied  an  ice  sheet  (at  C  band) 
that  was  sufficiently  smooth  that  volume  scatter  from  a 
6.5  cm  snow  cover  dominated  the  return  at  larger 
incidence  angles.  At  incidence  angles  greater  than  25° 
an  increase  of  up  to  7  dB  in  o°  was  observed. 


Figure  6.10:  Effect  of  6.5cm  cover  on  backscatter 

cross-section  of  smooth  saline  Ice  at  X 
band,  as  determined  from  data  obtained 
during  the  CRREL'87  experiment,  for  (a) 
VV  polarization  and  (b)  HH  polarization. 
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Figure  6.11:  Effect  of  7cm  snow  cover  on  backscatter 
cross-section  of  smooth  saline  Ice  at  C 
band,  as  determined  from  data  obtained 
during  the  CRUEL' 89  experiment  (see  data 
points  on  Figs.  6.7  and  6.13). 


Surprisingly,  a  10  dB  reduction  at  13°  incidence  was 
also  noticed  (See  Fig.  6.11).  The  layer  model  of  Ulaby 
et  al  [1982]  is  used  to  illustrate  the  snow  volume 
scatter  effect  at  larger  incidence  angles. 

This  layer  model  has  the  form: 

«°-Tlt«ds  +  °fce  /L>  1  +  °  f»ow  4u.r(  (6.12) 

where  T  =  transmission  coefficient  of  the  air-snow 
Interface 

L  =  exp(k  dsec(S) ) 

ofi  =  see  equation  (6.11) 


142 


40-1 - - - i - 1 - 1 - 1 - 1 - - 

10  IS  20  25  30  35  40  45 

Inc  Anglo  (dog) 


*  1.8  cm 

snow  dcpih  =  7  cm 
ffjc*  =  0.024  cm 
J,ct  =  1.2  cm 


—  Fit  to  sno» 
covorsd  ica 
data 

-  Fit  to  bora 
saline  ica 
data 

—  Predicted  ■ 
0.05  g/cu  ca 
snow 

—  Predicted  ■ 
0. 1  g/cu  ce 
enoe 


50 


Figure  6.12:  Backscatter  cross-section  of  snow  covered 
smooth  saline  ice  surface  at  C  band,  as 
determined  from  data  obtained  during  the 
CRREL'89  experiment,  compared  with 
predictions  of  the  layer  model  of  Ulaby 
et  al.  (19821 . 


r  ,  =  snow  surface  backscatter  cross- 

section 

Figure  6.12  shows  predicted  tr°  for  the  combined 
surface  using  pr=0.05  g/cm3  and  0.1  g/cra3  (measured 
value  of  pt  was  0.06  g/cmJ ) .  Notice  that  the  return  is 
well  predicted  for  p s  =  0.1  g/cmJ  except  at  13° 
incidence.  Figure  6.13  illustrates  the  relative 
contributions  to  a°  by  the  snow  surface,  snow  volume, 
and  ice  for  i»s=0.1  g.cm  .  Notice  that  the  measurement 
at  13°  is  the  same  as  that  predicted  for  snow 
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Figure  6.13:  Backscatter  cross-section  of  snow  covered 
smooth  saline  ice  surface  at  C  band,  as 
determined  from  data  obtained  during  the 
CRREL'89  experiment,  and  predicted 
relative  contributions  of  snow  surface, 
snow  volume,  and  ice. 


sur f ace+volume  scatter  only. 

The  reduction  at  13°  incidence  cannot  be  explained 
by  snow  surface  roughness.  Scattering  was  weak  due  to 
the  low  snow  dielectric  constant  ( <=  1.2).  Two 
possible  reasons  for  this  reduction  have  been 
considered,  neither  of  which  can  be  verified.  It  is 
possible  that  we  altered  the  properties  of  the  ice 
surface  in  the  process  of  removing  the  snow  (we  studied 
the  snow  covered  surface  prior  to  the  bare  ice 
surface).  However,  we  found  no  evidence  that  the  snow 


cover  modified  the  ice  surface  or  that  ve  altered  the 

ice  surface  in  the  process  of  shoveling  and  sweeping 
snow;  the  ice  surface  was  dry  and  naturally  smooth 
after  the  snow  was  removed  (the  conditions  at  the  time 
of  the  measurement  were  ice  salinity  of  3o/oo,  dry 
snow,  and  cold  temperatures  of  «-10°C) .  Further,  we 
had  little  difficulty  removing  the  snow  all  the  way 
down  to  the  ice  surface. 

Another  possible  explanation  is  that  the  snow  may 
have  been  acting  as  an  impedance  matching  layer. 
Normally,  a  thin  layer  of  snow  settles  to  a  nearly 
uniform  density  within  a  period  of  a  few  hours  after  it 
falls  [Gow,  1989],  We  studied  the  snow  covered  ice 
sheet  less  than  2  hours  after  the  snow  fell.  If  only 
partial  settling  had  occurred,  it  could  have  been  that 
the  density  of  the  snow  varied  as  a  function  of  depth, 
thus  forming  a  layer  with  a  tapered  impedance  (with 
impedances  varying  from  about  Z;^esp»te  /vTTT  to 
Zf jpjce/Vl.4-  Such  a  layer  could  have  caused  partial 
matching  over  the  wide  r.f.  bandwidth  of  our  radar 
(Tapered  impedances  are  commonly  used  to  match  circuits 
of  different  impedance  over  vide  bandwidths;  e.g., 
(Izadian,  1988]  and  [Ahmed,  1981]).  This  set  of 
circumstances  is  not  typical,  however  (Gow,  1989]. 


6.2.3  Rough  saline  ice 

During  the  CRREL'88  experiment  backscatter  from 
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Figure  6.14:  Comparison  of  backscatter  cross-sections 
of  smooth,  rough,  and  desalinated  ice  at 
C  band,  as  determined  from  data  obtained 
during  the  CRREL'88  experiment  (see  data 
points  on  Figs.  6.15  and  6.16). 


smooth  and  rough  ice  surfaces  va3  compared.  After 
measurements  of  the  smooth  ice  sheet  were  performed, 
1x2  cm  ice  prisms  were  dumped  on  the  ice  surface  and 
more  measurements  were  done.  The  roughness  caused  an 
increase  in  the  incoherent  return  at  larger  incidence 
angles,  as  would  be  expected  (See  Fig.  6.14). 

The  criteria  for  neither  the  small  perturbation 
nor  the  Kirchhoff  surface  scattering  theories  is 
satisfied  in  thi3  case,  i.e.,  the  rms  height  is  too 
large  for  the  former  (kos^).6  and  v'To/is# .  7 ) ,  and  the 
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Figure  6.15:  Backscatter  cross-section  of  rough  saline 
ice  sheet  at  C  band,  as  determined  from 
data  obtained  during  the  CRRBL'88 
experiment,  compared  with  predictions  of 
the  Small  Perturbation  Method. 


radius  of  curvature  is  too  small  for  the  latter. 
Nevertheless  the  Small  Perturbation  Method  predicts  the 
backscatter  measured,  over  the  entire  angular  range 
(See  Fig.  6.15)  . 

6.2.4  Desalinated  ice 

Backscatter  measurements  of  desalinated  ice  were 
made  at  C  band,  during  the  1988  experiments,  and  at  C 
and  Ku  band  during  1989.  During  1988  the  results  shown 
in  Fig.  6.14  were  obtained.  Notice  that  o°  for 
desalinated  ice  decreased  about  as  quickly  as  a 
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Figure  6.16:  Backscatter  cross-section  of  desalinated 
ice  sheet  at  C  band,  as  determined  from 
data  obtained  during  the  CRREL'88 
experiment,  compared  with  predictions  of 
the  Small  Perturbation  Method. 


function  of  incidence  angle  as  did  a°  for  smooth  saline 

•  ice,  up  to  30°  incidence,  implying  that  surface  scatter 

was  dominant.  The  flat  portion  at  incidence  angles 
greater  than  30°  may  have  been  due  to  volume  scatter; 

^  it  is  difficult  to  be  certain  because  radar  system 

effects  were  considerable  at  this  low  /  level. 

Notice  from  Fig.  6.16  that  the  Small  Perturbation 
Method  predicts  the  backscatter  measured,  over  the 

•  entire  angular  range. 

A  simple  experiment  was  done  on  this  sheet  to 
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study  the  volume  scattering  component.  Backscatter  was 

measured  as  a  function  of  ice  surface  temperature  at  a 
fixed  incidence  angle.  We  found  that  C-band  varied 
almost  Imperceptibly  as  free  water  formed  at  the  ice 
surface  (Fig.  6.17a)  confirming  that  the  volume  scatter 
contribution  was  very  small.  In  contrast  we  observed 
that  a°  at  X  band  decreased  markedly  as  the  ice  surface 
began  to  melt  (Fig.  6.17b)  indicating  that  volume 
scatter  was  significant  at  X  band. 

During  1989  we  measured  backscatter  from 
desalinated  ice  blocks  (See  Chpt .  3).  The  results  of 
the  backscatter  measurements  at  C  and  Ku  bands  are 
shown  in  Fig.  6.18.  Notice  that  o°(0)  for  Ku  band  is 
approximately  independent  of  0,  characteristic  of 
volume  scatter  (Kim,  1984a].  On  the  other  hand  crQ(6) 
at  C  band  falls  off  more  than  10  dB  from  15°  incidence 

O 

to  25  incidence,  as  is  characteristic  of  a  return  with 
a  large  surface  scatter  component. 

During  the  following  analysis  we  assume,  for 
simplicity,  that  the  ice  blocks  contained  spherical  air 
bubbles  of  radius  0.45  mm,  that  the  volume  fraction  of 
air  in  the  ice  was  20%,  and  the  dielectric  constant  of 
the  ice  was  uniform  at  2.5-j0. 01;  these  were 
characteristics  of  the  top  layer  of  the  ice.  One-way 
transmission  losses  of  the  ice  at  Ku  and  C  bands  were 
10  dB  at  Ku  band  and  5  dB  at  C  band.  Predicted  one-way 
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Figure  6.18:  Backscatter  cross-section  of  desalinated 
ice  at  (a)  C  band  and  (b)  Ku  band,  as 
determined  from  data  obtained  during  the 
CRREL'89  experiment,  compared  vith 
combined  predictions  of  Small 
Perturbation  Method  surface  and  Effective 
Medium  volume  scattering  theories. 


attenuation  due  to  absorption  for  this  ice  sheet,  based 
on  the  estimated  dielectric  profile  is  shown  in  Fig. 
3.21  and  is  3.5  dB  at  Ku  band  and  1.5  dB  at  C  band. 
Thus,  scattering  losses  are  assumed  to  be  6.5  dB  at  Ku 
band  and  3.5  dB  at  C  band.  Because  the  air  bubble 
sizes  were  larger  and  the  number  of  them  more  dense  in 
the  upper  layer,  it  is  assumed  that  they  were 
responsible  for  the  greatest  scattering  loss  and  that 
we  are  justified  in  making  the  above  assumptions. 

As  a  first  approximation  we  assume  that  the 
backscatter  was  due  to  uncorrelated  Rayleigh 
scattering.  The  volume  scatter  coefficient  in  this 
case  is 

av  =  64TT5/*4ice  IK1  iNr4*  [Ulabyetal,  198  ]  (6.13) 

where  *ice  =  wavelength  in  ice 

K  =  (  «c  -1 )  /  (  +  2  )  =  0.33  (where  €c=l/«;te  ) 

N  =  number  of  air  bubbles/cu  meter 
r  =  radius  of  air  bubbles 
At  C  band  av  a;  0.016  and  at  Ku  band  ov  a;  0.647 
Assuming  that  the  extinction  coefficient  is  uniform 
throughout  the  ice  sheet  and  is  of  magnitude  to  account 
for  the  losses  measured,  ke  =  -7.2  for  C  band  and  -14.4 
for  Ku  band. 

« 

This  results  in  o°  a:  ov  Jexp  ( -2kt  z  )dz  (6.14) 

0 

or  o®  aj  -16.5  dB  at  Ku  band  and  -29.5  dB  at  C  band. 
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These  results  are  considerably  lower  than  the  data 
(See  Fig.  6.18  for  data). 

Consider  predictions  of  the  Effective  Medium 
volume  scattering  theory.  This  theory  considers  both 
correlated  and  uncorrelated  scattering.  (The  theory 
that  was  used  differs  slightly  from  the  one  presented 
in  Chapter  2  -  we  used  the  one  considered  here  because 
tested  software  was  already  available.  This  theory 
differs  from  the  one  presented  in  Chapter  2  in  that 
this  theory  only  considers  the  loss  due  to  scattering 
that  is  accounted  for  in  estimates  of  the  dielectric 
constant,  and  correlation  between  scatterers  is 
considered  by  means  of  a  correlation  function  rather 
than  a  two-particle  distribution  function.  Details 
can  be  found  in  Fung  [1989]  .  )  The  predictions  for 
the  volume  scatter  component  of  the  desalinated  ice 
blocks,  using  volume  fraction  of  air  =  20%,  air  bubble 
radius  =  0.45  mm,  correlation  length  =1.5  mm,  and  e" 

=  0.01,  and  Small  Perturbation  predictions  for  the 
surface  scatter  component,  are  shown  in  Figs.  6.18. 
Notice  that  the  prediction  for  Ku  band  falls  within  the 
measurement  error  bars.  On  the  other  hand  the  C  band 
prediction  is  7  or  8  dB  higher  than  the  measured  data. 
We  are  uncertain  of  the  reason  for  the  large  prediction 
at  C  band. 
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To  determine  a  possible  explanation  for  this  ve 
varied  each  input  parameter  around  the  measured  value 
to  study  its  effect.  Notice  from  Fig.  6.19  that  a 
decrease  in  volume  fraction  causes  the  prediction  at 
both  frequencies  to  decrease.  The  same  is  true  of  an 
increase  in  «"  (Fig.  6.20).  Notice  in  Fig.  6.21, 
however,  that  a  decrease  in  the  correlation  length 
caused  a  large  decrease  in  C  band  predictions  but  only 
a  small  decrease  in  the  Ku  band  predictions. 

The  results  of  these  comparisons  of  measurements 
of  desalinated  ice  blocks  and  theoretical  predictions 
are  encouraging;  but  several  things  prevent  us  from 
making  stronger  statements  concerning  quantitative 
comparisons.  First,  the  bubble-structure  data  were 
obtained  from  one  sample  of  the  desalinated  ice  blocks, 
leaving  questions  as  to  the  statistical  uncertainty  of 
the  bubble  structure  -  Figure  6.21  shows  that  such 
uncertainty  can  be  of  considerable  importance.  Second, 
the  conditions  under  which  the  experiments  were 
performed  were  iess  than  ideal.  For  example,  the 
antenna  pointing  angle  was  difficult  to  control  (and 
monitor)  for  this  experiment  because  the  gantry  was  on 
a  very  unlevel  surface,  leaving  questions  as  to  the 
uncertainty  in  surface  contributions  at  lower  incidence 
angles;  and  It  was  not  possible  to  eliminate  or 
quantify  the  effects  of  backscatter  from  block  edges. 
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Figure  6.19:  Backscatter  cross-section  of  desalinated 
Ice  at  (a)  C  band  and  (b)  Ku  band,  as 
determined  from  data  obtained  during  the 
CRREL'89  experiment,  compared  with 
predictions  of  the  Effective  Medium 
theory  at  three  different  values  of 
volume  fraction  of  air. 
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Figure  6.20:  Backscatter  cross-sect  Ion  of  desalinated 
Ice  at  (a)  C  band  and  (b)  Ku  band,  as 
determined  from  data  obtained  during  the 
CRREL'89  experiment,  compared  with 
predictions  of  the  Effective  Medium 
theory  for  three  different  values  of  the 
imaginary  component  of  the  dielectric 
constant . 
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Figure  6.21:  Backscatter  cross-section  of  desalinated 
ice  at  (a)  C  band  and  (b)  Ku  band,  as 
determined  from  data  obtained  during  the 
CRREL'89  experiment,  compared  with 
predictions  of  the  Effective  Medium 
theory  for  three  different  values  of 
correlation  length  of  spacing  between  air 
bubbles . 
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which  could  have  been  considerable  since  vide-beamvidth 


antennae  (horns)  were  used,  resulting  in  uncertainty  as 
to  the  relative  contributions  of  both  surface  and 
volume . 
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The  purpose  of  this  research  was  to  determine  the 
sources  of  scattering  of  artificial  sea  ice,  to  compare 
backscatter  measurements  semi-quant itatively  with 
theoretical  predictions,  and  to  develop  inexpensive 
polarimetric  radars  for  sea  ice  backscatter  studies. 

In  the  investigation  of  scattering  from  artificial 
sea  ice  we  came  to  the  following  conclusions: 

1.  The  volume  scatter  contribution  to  backscatter  from 
smooth  (rms  height  <  0.05  cm)  saline  ice  at  C  band 
cannot  be  neglected.  We  were  unable  to  quantify  the 
relative  surface  and  volume  contributions,  but  our 
results  indicate  that  the  volume  contribution  is 
comparable  to  the  surface  contribution. 

2.  A  thin  dry  snow  layer  greatly  altered  the  measured 
backscatter  cross-section  of  smooth  (rms  height  < 
0.05  cm)  saline  ice  at  C  band  and  incidence  angles 
greater  than  about  25° .  The  large  reduction  in  o>° 
near  nadir,  presented  in  Chapter  6,  was  due  to  a  set 
of  circumstances  atypical  of  conditions  in  the 
Arctic  or  Antarctic.  A  thin  dry  snow  layer  did  not 
noticeably  influence  measured  backscatter  f ? om 
slightly  rough  (rms  height  >  0.1  cm)  saline  ice  at  X 
band  when  the  surface  was  sufficiently 

rough  with  respect  to  a  wavelength  that  surface 
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3. 


4. 


scatter  dominated  over  volume  scatter  from  the  snow. 
One  can  surmise  that  snow  would  dominate  X-band 
scattering  if  the  underlying  ice  were  sufficiently 
smooth,  but  no  such  experiment  was  performed. 
Backscatter  from  rough  saline  ice  at  C  band  is 
dominated  by  the  contribution  from  the  surface.  For 
this  rough  surface  the  Small  Perturbation  Method 
correctly  predicted  the  measured  data  over  incidence 
angles  of  0  to  60  degrees. 

Backscatter  from  slightly  rough  (rms  height  >  0.15 
cm)  desalinated  ice  at  C  band  and  incidence  angles 
less  than  about  30  degrees  is  dominated  by  surface 
scatter.  Volume  scatter  is  important  at  X  band  and 


is  dominant  at  Ku  band. 

The  CRREL  facility  is  useful,  but  a  number  of 
factors  limited  the  quality  of  our  measurements.  These 
need  to  be  improved;  e.g., 

(1)  environmental  conditions  were  often  not  well 
controlled  because  of  the  poor  condition  of  the  roll¬ 
back  tent; 

(2)  a  common  gantry  restricted  the  experiment  plans  of 
Individual  investigators,  and  the  gantry  could  not 
easily  be  moved  by  a  single  person; 

(3)  insufficient  manpower  was  available  to  thoroughly 
characterize  the  ice  sheets;  and 

(4)  available  technology  was  not  being  used  to 
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characterize  the  ice.  For  example,  surface  roughness 
measurements  were  made  using  cumbersome  and  time- 
consuming  contact  methods  that  limited  the  number  of 
roughness  samples  that  could  reasonably  be  taken, 
rather  than  noncontact  laser  profilometer  or  stereo 
pair  photographic  techniques;  and  dielectric  constants 
were  estimated  from  empirical  equations  rather  than 
measured . 

The  uncertainty  in  our  backscatter  data  could  have 
been  reduced  if  we  would  have  been  able  to  carry  out 
our  own  experiment  plans  because  in  many  cases  we  would 
have  taken  more  measurements  than  we  were  otherwise 
permitted  to,  and  the  uncertainty  in  surface 
characterization  could  have  been  reduced  if  available 
technology  and  more  manpower  had  been  utilized  to 
thoroughly  characterize  the  ice  sheets. 

Finally,  we  conclude  that  inexpensive,  radars,  such 
as  those  presented  in  Chapter  4,  can  provide  adequate 
measurements  of  backscatter  cross-section  and 
scattering  matrices  of  the  surfaces  reported  in  this 
paper . 

7.2  Recommendations 

Future  measurements  should  emphasize  thorough 
characterization  of  the  ice  surfaces  so  that  more 


precise  comparisons  of  measurement  and  theory  can  be 
done.  Sufficient  manpower  should  be  provided  to 


continuously  monitor  ice  structural  parameters;  and 
available  technology  should  be  utilized  to  measure 
these  parameters  frequently  and  accurately.  Current 
technology  is  being  used  to  study  the  internal 
characteristics  of  the  ice.  It  should  also  be  used  to 
monitor  surface  roughness,  dielectric  properties,  etc. 

Measurements  of  very  smooth  surfaces  should  be 
conducted  with  lov-sidelobe  narrow-beam  antennas  to 
improve  sensitivity  and  to  avoid  problems  introduced  by 
vide-beam  antennas. 

The  phase  calibration  technique  that  was  developed 
is  useful  in  the  laboratory.  To  avoid  the  need  for 
frequent  phase  calibration,  stable  oscillators  should 
be  used  for  up-converting  and  down-converting  in  step 
frequency  radars,  and  VV  and  HH  polarization  channel 
system  paths  should  be  phase-matched,  to  the  extent 
possible,  for  the  FM-CW  and  step  frequency  radar 
systems . 
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